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In  support  of  the  Defense  Nuclear  Agency  (DNA)  Silo  Test  Program,  the  planar  steady  detona¬ 
tion  processes  were  characterized  and  initial  JWL  (Jones-Wilkins-Lee)  equation  of  state 
(EOS)  coefficients  were  developed  for  two  chemical  high  explosives  (HEs),  NA/NP  and  NPN. 
NA/NP  is  a  liquid  explosive  composed  of  weak  (Baum€  42)  aqueous  nitric  acid  (NA)  and  the 
solvent  nitropropane  (NP).  The  NA/NP  studied  here  was  3.158  parts  NA  to  1  part  NP,  by 
weight,  and  its  density  was  1.25  Mg/m^  at  29®G.  NPN  (nitropropane  nitrate)  is  a  damp  gran¬ 
ulated  blasting  agent  composed  of  ammonium  nitrate  (AN),  NP,  and  methyl  alcohol  (MA).  Our 
study  material  was  0.848  AN,  0.075  NP,  0.075  MA,  and  0.002  methocel  (an  added  antimigratory 
agent),  by  weight,  and  had  a  density  of  1.01  Mg/m’.  These  unusual  explosives  are  being 
considered  by  DNA  for  use  in  various  nuclear  weapons  effects  simulators  because  of  their 
potential  for  producing  desired  pulse  shapes  at  competitive  costs. 

SRI  large-scale  multiple  Lagrange  particle  velocity  gage  experiments,  in  combination  with 
Lagrange  analyses,  quantified  the  early  stages  of  the  steady  planar  detonation  process  in 
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19.  ABSTRACT  (Continued) 

each  material.  Thermodynamic  equilibrium  calculations,  performed  with  the  SRI  TIGER  code, 
extended  the  adiabatic  release  paths  to  pressures  approaching  atmospheric. 

NA/NP  exhibited  ideal  detonation  with  a  pressure-relative  volume  (p-V)  adiabatic  release 
path  slightly  above  that  of  nitromethane  (NM)  and  a  resolved  reaction  zone  of  surprisingly 
long  duration,  0.2-0. 6  ys.  A  single  experiment  indicated  that  the  shock-to-detonation 
transition  (SDT)  may  occur  over  unusually  long  buildup  distances. 

The  NPN  detonation  process  was  nonideal  in  that  not  all  the  AN  reacted  to  support  the 
steady  detonation.  The  p-V  release  adiabat  was  very  close  to  that  for  ANFO  over  a  large 
relative  volume  range,  but  lay  somewhat  above  it  at  and  near  the  detonation  state.  (At 
higher  initial  densities,  the  NPN  performance  advantages  near  the  detonation  state  would  be 
greater.)  Two  possible  two-step  AN  reactions  were  postulated  to  model  the  nonideal  detona¬ 
tion  process  in  NPN. 

A  new  JWL  fitting  procedure  was  developed  based  on  the  detonation  state  release  isentrope 
and  three  hydro-thermodynamic  constraints.  The  procedure  is  intended  to  provide  JWL  EOS 
parameters  suitable  for  use  in  DNA  hydrocode  calculations  to  estimate  high  explosive 
performance  in  various  nuclear  weapons  effects  simulators.  The  procedure  was  used  to  deter 
mine  JWL  EOS  parameters  for  each  material. 

Five  recommendations,  three  strong,  are  made  to  DNA  in  the  areas  of  HE  use  and  characteriza 
tion: 

(1)  NA/NP  is  a  very  promising  liquid  HE;  we  strongly  recommend  considering  its  use  in  NWE 
simulators,  as  a  low  cost  alternative  to  NM. 

(2)  Although  NPN  compacts  irreversibly,  especially  at  low  densities,  we  recommend  its 
consideration  in  NWE  applications  calling  for  an  HE  similar  to  ANFO  but  with  the  ability  to 
do  somewhat  more  work  at  high  pressures  and/or  with  more  rapid  kinetics. 

(3)  We  strongly  recommend  Lagrange  gage  SDT  experiments  to  characterize  the  buildup  and 
propagation  of  detonation  in  NA/NP. 

(4)  We  recommend  further  research  to  confirm  the  presumed  charge  size  independence  of  the 
present  NPN  experimental  results. 

(5)  We  strongly  recommend  further  experimental  and  computational  efforts  to  improve  the 
understanding  of,  and  ability  to  model,  nonideal  detonation  processes  in  AN-rich  HEs  such 
as  NPN  and  ANFO. 
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Defense  Nuclear  Agency  (DNA)  for  possible  use  in  various  nuclear  weapon 
effects  (NWE)  simulators  in  the  Silo  Test  Program  and  other  programs. 
These  novel  HEs  may  be  able  to  produce  desired  pulse  shapes  at  competi¬ 
tive  cost  in  such  applications.  In  the  present  research,  SRI 
International  Investigated  NA/NP  and  NPN  to  (1)  experimentally  and  com¬ 
putationally  characterize  the  planar  steady  detonation  process  in  each 
HE  at  large  charge  diameters  and  (2)  develop  initial  JWL  (Jones-Wilkins- 
Lee)  equation  of  state  (EOS)  parameters  suitable  for  DNA  hydrocode 
calculations  predicting  the  performance  of  each  HE  in  various  NWE 
simulators . 

METHODS 

SRI  large-scale  multiple  Lagrange  particle  velocity  gage  experi¬ 
ments,  in  combination  with  Lagrange  analyses,  were  used  to  quantify  the 
early  stages  of  the  steady  planar  detonation  processes  in  each  material. 
In  the  experimental  technique,  thin  metallic  foil  gages  were  embedded  at 
various  depths  within  the  explosive  to  provide  the  velocity  histories  of 
particles  (Lagrange  histories)  during  the  shock  compression  and  early 
expansion  phases  of  the  detonation  process. 

Thermodynamic  equilibrium  calculations  were  used  to  complete  the 
quantification  of  the  detonation  process  by  extending  the  adiabatic 
expansion  path  determined  experimentally  to  pressures  approaching  atmos¬ 
pheric.  These  calculations  assume  chemical,  mechanical,  and  thermal 
equilibrium  among  the  detonation  products  and  determine  the  composition 
and  state  variables  for  proscribed  thermodynamic  processes  modeling  the 
adiabatic  expansion.  The  calculations  were  performed  with  the  SRI  TIGER 
code  and  used  the  BKW  and  JCZ3  EOSs  for  individual  detonation  product 
gases. 
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A  new  fitting  procedure  was  developed  to  determine  JWL  EOS  parame¬ 
ters  from  the  adiabatic  release  path  through  the  CJ  (Chapman-Jouguet) 
detonation  state.  The  procedure  provides  unique  values  for  the  JWL 
parameters  by  minimizing  the  RMS  deviation  between  the  CJ  release 
adlabat  and  the  JWL  CJ  Isentrope  subject  to  three  hydro-thermodynamic 
constraints.  The  procedure  Is  quite  general  and  Is  believed  to  have 
applications  beyond  the  present  program.  Initial  JWL  parameters  were 
developed  for  each  material  by  applying  this  procedure  to  the  release 
adlabat  determined  from  the  Lagrange  gage  experiments  and  extended  by 
the  TIGER  calculations. 

NA/NP  CHARACTERIZATION 

NA/NP  Is  a  liquid  explosive  composed  of  Baume  42  nitric  acid  (NA), 
a  weak  aqueous  nitric  acid  solution  containing  67.18  wt%  pure  nitric 
acid,  and  nltropropane  (NP),  a  solvent.  The  NA/NP  composition  studied 
here  was  3.158  parts  NA  to  1  part  NP,  by  weight,  and  had  a  density  of 
1.25  Mg/m^  at  29®C.  This  composition  Is  oxygen  balanced  to  CO2. 

Despite  Its  acid  and  solvent  components,  we  found  NA/NP  to  be  relatively 
easy  to  handle. 

Two  successful  Lagrange  gage  detonation  experiments  were  performed 
In  NA/NP.  The  second  had  more  precise  Instrumentation  settings  and  pro¬ 
duced  high  quality  reproducible  particle  velocity  data  that  were 
Lagrange  analyzed  to  determine  the  compression  and  initial  adiabatic 
release  states  of  a  particle  In  the  steady  detonation  process.  NA/NP 
exhibited  Ideal  detonation  conforming  to  the  ZND  (Zeldovlch,  Von 
Neumann,  Doering)  model  with  a  clearly  resolved  reaction  zone  of  sur¬ 
prisingly  long  duration,  0.2-0. 6  ps.  At  the  CJ  state  the  pressure  was 
13.3  GPa,  the  detonation  velocity  was  6.48  km/s,  the  particle  velocity 
was  1.64  km/s,  the  specific  volume  was  0.5986  m^/Mg,  and  the  adiabatic 
exponent  was  2.99. 

A  single  Lagrange  gage  experiment  with  a  reduced  amplitude  Input 
shock  was  also  performed  to  examine  the  shock-to-detonatlon  transition 
(SDT)  In  NA/NP.  The  resulting  particle  velocity  histories,  although  not 


reproducible  across  the  target,  indicate  the  possibility  of  an  unusual 
SDT  phenomenology  with  long  buildup  distances.  If  confirmed  by  further 
experiments,  this  effect  should  be  accounted  for  In  NWE  simulator 
designs  using  NA/NP. 

To  extend  the  CJ  adiabatic  release  path  from  the  second  Lagrange 
gage  NA/NP  detonation  experiment  to  large  volumes,  we  assumed  the  expan¬ 
sion  to  be  Isentroplc  and  used  TIGER  with  both  the  BKW  and  the  JCZ3  EOSs 
to  estimate  the  Isentrope.  The  BKW  calculations  terminated  prematurely 
at  about  1  GPa,  whereas  the  JCZ3  calculations  extended  to  about  2  MPa. 
Therefore,  the  JCZ3  calculations  were  selected  to  characterize  the  NA/NP 
detonation  product  expansion  to  large  volumes.  The  TIGER/JCZ3  Isentrope 
through  the  experimental  Cl  state  agreed  very  closely  with  the  release 
adlabat  from  the  Lagrange  gage  experiment  in  the  regions  of  overlap. 

NA/NP  JWL  EOS  coefficients  were  calculated  by  applying  the  new 
fitting  method  to  the  CJ  state  determined  from  the  Lagrange  gage  detona¬ 
tion  experiment  and  the  release  path  through  It  from  the  TIGER/JCZ3 
calculations.  The  resulting  coefficients  are  as  follows: 

A  =  7,554  Mbar,  B  =  0.1962  Mbar,  C  -  0.012717  Mbar,  Rj  *  6.366, 

R2  =  1.810,  and  w  =  0.285.  These  coefficients  are  the  major  result  of 
the  present  NA/NP  characterization  effort.  They  are  judged  appropriate 
for  use  in  hydrocode  calculations  to  evaluate  NA/NP  performance  in  a 
variety  of  NWE  simulator  applications. 

The  NA/NP  JWL  CJ  release  isentrope  determined  here  is  close  to,  but 
slightly  above,  that  for  nitromethane  (NM)  determined  from  cylinder  test 
data.  Therefore,  NA/NP  performance  in  NWE  simulators  Is  expected  to  be 
very  similar  to  that  of  NM,  although  its  material  cost  Is  much  less.  We 
conclude  that  NA/NP  Is  a  very  promising  liquid  HE  offering  an  attractive 
low  cost  alternative  to  NM. 

NPN  CHARACTERIZATION 

NPN  (nltropropane  nitrate)  is  a  damp  granulated  blasting  agent 
composed  of  ammonium  nitrate  (AN),  NP,  and  methyl  alcohol  (MA).  The 
NPN  studied  here  was,  by  weight,  0.848  ground  porous  AN  prills,  0.075 


NP,  0.075  MA,  and  0.002  methocel  (an  added  anti  ory  agent)  and  had 

a  density  of  1.01  Mg/m^.  This  composition  is  very  slightly  overfueled 
relative  to  CO2  balance.  The  NPN  was  selected  to  be  close  to  that 
examined  previously  for  DNA  In  an  alrblast  experiment.  We  found  that 
careful  handling  was  required  to  keep  NPN  In  a  uniform  state  because  It 
easily  undergoes  Irreversible  compaction  at  these  low  densities. 

Special  procedures  were  developed  to  uniformly  load  the  NPN  in  the 
Lagrange  gage  experiment  target  chambers  without  distorting  the  embedded 
particle  velocity  gages.  Two  NPN  multiple  Lagrange  particle  velocity 
gage  experiments  were  performed,  but  despite  the  precautions,  some  gage 
displacement  occurred  In  the  first  experiment.  Therefore,  only  data 
from  the  second  experiment.  In  which  slightly  modified  loading  proce¬ 
dures  were  used,  were  analyzed  quantitatively. 

For  the  valid  recording  duration  of  the  experiment  (about  8  ps), 
the  reactive  flow  In  the  NPN  consisted  of  a  uniform  particle  velocity 
pulse  propagating  at  a  velocity  of  5.08  km/s.  The  steady  particle 
velocity  history  profile  Is  Idealized  as  consisting  of  a  shock  compres¬ 
sion  to  a  peak  amplitude  of  about  2.4  km/s  In  0.2  ps,  a  decay  to  about 
0.6  of  the  peak  amplitude  In  about  2  ps,  and  a  very  slow  further  decay, 
or  a  constant  state,  for  about  6  ps.  The  steady  flow  and  slow  final 
decay  rate  observed  in  the  NPN  have  not  been  observed  in  the  other  HEs 
we  have  studied  with  this  technique  and  suggest  that  the  detonation  pro¬ 
cess  in  NPN  may  be  nonideal. 

Because  the  flow  observed  In  the  NPN  detonation  experiment  was 
steady,  the  Lagrange  analysis  could  be  performed  by  application  of  the 
Ranklne-Hugonlot  relations  for  shocks  and  steady  flows.  The  resulting 
pressure  history  has  the  same  shape  as  the  Idealized  particle  velocity 
history,  and  the  pressure-relative  volume  path  consists  of  compression 
and  expansion  along  the  Rayleigh  line. 

An  extensive  effort  was  required  to  simulate  and  extend  the  experi¬ 
mental  results  with  TIGER.  Initial  calculations  performed  with  both  the 
BKW  and  the  JCZ3  EOSs  revealed  that  the  detonation  process  was  indeed 
nonideal  In  that  not  all  the  AN  reacted  to  support  the  detonation.  BKW 


predicted  that  51.5  wt%  of  the  AN  did  not  support  the  detonation  and 
estimated  the  CJ  pressure  at  5.44  GPa.  JCZ3  predicted  that  only  14.5 
wt%  of  the  AN  was  Inert  and  estimated  the  CJ  pressure  at  6.48  GPa. 

Because  of  the  different  calculated  amounts  of  unreacted  AN  at  the 
CJ  state  and  because  a  BKW  procedure  had  been  developed  previously  for 
treating  nonideal  detonation  in  ANFO  (an  HE  similar  to  NPN),  we  postu¬ 
lated  different  secondary  AN  reaction  processes  to  calculate  adiabatic 
release  paths  to  large  volume  with  the  two  EOSs.  In  the  BKW  calcula¬ 
tions,  following  the  previous  treatment  for  ANFO,  we  assumed  that  the 
remaining  51.5  wt%  of  AN  was  consumed  in  an  irreversible  reaction  spread 
over  a  large  pressure-volume  range  starting  at  the  CJ  state  and  that  the 
subsequent  expansion  was  isentroplc.  However,  in  the  JCZ3  calculations, 
because  of  the  nearly  constant  state  observed  in  the  latter  part  of  the 
experimental  records,  we  assumed  that  the  small  remaining  amount  of  AN 
reacted  slowly  at  a  pressure  and  volume  near  that  of  the  CJ  state  and 
that  the  subsequent  expansion  was  isentroplc.  Because  only  the  JCZ3 
results  extended  to  large  enough  relative  volumes  to  permit  application 
of  the  new  JWL  fitting  procedures,  we  used  the  JCZ3  calculations  to 
represent  the  NPN  CJ  release  adiabat. 

Initial  JWL  EOS  parameters  were  calculated  for  NPN  by  applying  the 
new  fitting  procedures  to  the  CJ  release  adiabat  consisting  of  the 
Rayleigh  line  determined  from  the  Lagrange  gage  experiment  and  the  isen- 
trope  through  its  final  state  from  the  TIGER/JCZ3  calculations.  The 
resulting  coefficients  are  as  follow:  A  =  6.783  Mbar,  B  =  0.1824  Mbar, 

C  =  0.007212  Mbar,  Rj^  =  7.628,  R2  =  2.305,  and  w  *  0.345.  These  coeffi- 
.-ients  are  judged  appropriate  for  use  in  hydrocode  calculations  to 
evaluate  NPN  performance  in  a  variety  of  NWE  simulator  applications. 

The  NPN  JWL  release  isentrope  determined  in  this  work  is  very  close 
to  that  for  ANFO  from  cylinder  tests,  over  the  full  relative  volume 
range  examined,  but  lies  somewhat  above  it  at  and  near  the  detonation 
state.  At  higher  initial  densities,  the  NPN  performance  advantages  over 
ANFO  would  be  expected  to  be  greater.  The  enhanced  high  pressure  per¬ 
formance  of  NPN  may  outweigh  its  handling  and  cost  disadvantages  rela¬ 
tive  to  ANFO  in  some  DNA  applications. 


The  nonideal  AN  reaction  phenomenology  is  presumed  to  be  similar  in 
NPN  and  ANFO,  but  the  findings  of  this  study  indicate  that  this 
phenomenology  is  not  yet  defined  or  quantified*  Therefore,  the  ability 
to  accurately  predict  the  performance  of  either  material  in  specific 
applications  is  suspect*  This  is  an  important  conclusion  since  ANFO 
will  probably  be  the  HE  used  by  DNA  in  forthcoming  large  airblast  simu¬ 
lations* 

RECOMMENDATIONS 

The  results  of  the  NA/NP  and  NPN  characterization  efforts  lead  us 
to  make  five  recommendations  to  DNA  in  the  areas  of  HE  use  and  charac¬ 
terization*  Recommendations  (1),  (3),  and  (5)  are  recommended  strongly* 


(1)  We  strongly  recommend  the  use  of  NA/NP  in  NWE  simulation  appli¬ 
cations  calling  for  low  cost  liquid  explosives  similar  in 
performance  to  NM* 

(2)  We  recommend  that  the  use  of  NPN  be  considered  in  NWE  applica¬ 
tions  calling  for  an  HE  similar  to  ANFO  but  with  the  ability  to 
do  somewhat  more  work  at  high  pressures  and/or  with  more  rapid 
kinetics*  Its  relative  performance  advantages  over  ANFO,  and 
the  ease  of  keeping  it  uniform,  probably  both  Increase  at 
higher  Intlal  densities  than  studied  here* 

(3)  We  strongly  recommend  support  for  Lagrange  gage  shock-to-deto- 
natlon  transition  experiments  to  characterize  the  buildup  and 
propagation  of  detonation  in  NA/NP. 

(4)  We  recommend  further  research  to  confirm  the  presumed  charge 
size  Independence  of  the  present  NPN  experimental  results* 

(5)  We  strongly  recommend  support  for  further  experimental  and 
computational  efforts  to  improve  the  understanding  of  and 
ability  to  model,  nonideal  detonation  processes  in  AN-rich  HEs 
such  as  NPN  and  ANFO  * 
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SECTION  1 


INTRODUCTION 


The  Defense  Nuclear  Agency  (DNA)  Is  using  chemical  high  explosives 
(HEs)  to  simulate  various  nuclear  weapon  effects  (NWEs)  Including  direct 
and  alrblast  Induced  ground  shock  and  alrblast  effects  on  structures* 

In  such  simulations,  kilotons  of  HE  and/or  very  precise  blast  pulse 
shapes,  amplitudes,  and  durations  may  be  required.  These  simulations 
have  generated  the  need  for  new  HEs  offering  lower  cost.  Improved  per¬ 
formance,  or  different  properties  than  conventional  HEs  previously 
developed  for  military  or  mining  purposes. 


In  this  work  we  examined  two  new  HEs  of  Interest  to  DNA:  a  liquid 
explosive  here  called  NA/NP  and  a  blasting  agent  called  NPN.  The  over¬ 
all  objective  was  to  generate  initial  JWL  (Jones-Wllklns-Lee)  equatlon- 
of-state  parameters  for  each  HE  to  allow  hydrocode  calculations  for 
evaluating  the  performance  of  each  material  In  NWE  simulators. 


NA/NP  Is  mixture  of  aqueous  nitric  acid  (NA)  and  the  solvent  nitro- 
propane  (NP).  NA/NP  offers  performance  similar  to  nitromethane  (NM)  at 
a  much  lower  material  cost.  NPN  (nltropropane  nitrate)  Is  composed  of 
ammonium  nitrate  (AN),  NP,  and  methyl  alcohol  (MA).  NPN  is  similar  to 
the  commercial  blasting  agent  ANFO  (ammonium  nitrate  and  fuel  oil)  used 
in  the  mining  Industry,  but  has  a  higher  detonation  velocity  and  pres¬ 
sure  and  a  somewhat  greater  cost.  Both  materials  are  being  studied  for 
the  DNA  Silo  Test  Program,  but  have  potential  for  many  other  DNA  appli¬ 


cations  as  well. 


The  present  work  included  both  experimental  and  theoretical  charac¬ 
terization  techniques.  The  SRI  large-scale  Lagrange  particle  velocity 
gage  method  was  used  to  measure  particle  velocity  histories  within  each 
HE  as  it  underwent  steady  detonation.  Lagrange  analyses  of  these  data 
provided  histories  of  the  other  flow  variables  including  propagation 
velocity,  pressure,  specific  volume,  and  specific  Internal  energy  during 


the  chemical  reaction  and  the  initial  adiabatic  expansion  of  the  detona¬ 
tion  products  for  each  material*  Thermodynamic  equilibrium  calculations 
performed  with  the  SRI  TIGER  code  modeled  the  Lagrange  gage  and  analysis 
results  and  extended  the  adiabatic  product  expansion  data  to  large 
^  volumes.  Initial  JWL  coefficients  suitable  for  hydrocode  calculations 

were  developed  for  each  material  from  the  results  of  the  Lagrange  gage 
'  experiments  and  the  TIGER  calculations. 

Section  2  of  this  report  describes  the  techniques  we  used,  includ¬ 
ing  a  valuable  new  procedure  for  obtaining  JWL  coefficients  consistent 
J  with  a  specified  detonation  state  and  the  release  isentrope  passing 

I  through  it.  Sections  3  and  4  present  the  results  of  the  NA/NP  and  NPN 

characterization  efforts,  respectively.  Section  5  gives  our  recommenda¬ 
tions  for  the  use  and  further  characterization  of  NA/NP  and  NPN. 


SECTION  2 


METHODS 

Three  methods  were  used  to  quantify  the  detonation  state  and  the 
subsequent  adiabatic  expansion  states  generated  in  HE  reaction  products 
by  steady  planar  detonations:  the  SRI  large-scale  Lagrange  particle 
velocity  gage  technique,  GUINSY  uniaxial  flow  Lagrange  analysis  calcu¬ 
lations,  and  TIGER  thermodynamic  equilibrium  calculations.  These 
methods  are  reviewed  below.  A  new  technique  developed  to  calculate 
unique  JWL  equatlon-of-state  (EOS)  parameters  from  the  detonation  state 
and  the  release  isentrope  centered  on  it  is  described  here  for  the  first 
time . 

2.1  SRI  LARGE-SCALE  LAGRANGE  PARTICLE  VELOCITY  GAGE  TECHNIQUE 

The  objective  of  this  technique  is  to  measure  histories  along 
particle  paths  (Lagrange  histories)  of  the  particle  velocity  at  several 
depths  within  an  explosive  undergoing  steady  planar  detonation.  Here 
and  throughout  this  report  we  use  as  the  basis  for  discussion  the  well- 
known  ZND  (Zeldovich-VonNeumann-Doering)  model  for  ideal  steady  planar 
detonation  with  an  explicit  reaction  zone.  In  this  model  the  detonation 
process  experienced  by  a  particle  is  idealized  as  consisting  of  (1) 
nonreactlve  shock  compression,  resulting  from  a  shock  propagating  at  the 
detonation  velocity  D,  to  a  state  on  the  Hugonlot  of  the  unreacted 
explosive  ("spike"  state),  (2)  chemical  reaction  and  expansion,  produced 
by  a  steady  rarefaction  (steady-state  reaction  zone)  also  propagating  at 
the  detonation  velocity  D,  terminating  at  a  state  on  the  detonation 
products  Hugonlot  (CJ  [Chapman-Jouguet]  detonation  state),  and  (3) 
nonreactlve  isentropic  expansion  to  ambient  pressure,  resulting  from  a 
dispersive  rarefaction  (Taylor  wave)  in  the  detonation  products. 

The  representation  of  the  ideal  ZND  detonation  process  in  the 
pressure-specific  volume  (p-v)  plane  and  the  corresponding  Lagrange 
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particle  velocity  (u)  histories  are  shown  schematically  In  Figure  1. 

The  objective  of  the  present  Lagrange  particle  velocity  gage  experiments 
Is  to  measure  the  actual  particle  velocity  histories  analogous  to  those 
In  Figure  1  so  that  the  real  detonation  process  and  associated  parame¬ 


ters  can  be  determined. 


HE  Hugoniot 
Products  Hugoniot 
Peak  Compression  State 
CJ  State 

1  Nonreactive  Shock 
Compression  Path 

2  Reaction  States 

3  Products  CJ  Release 
Isentrope 


SPECIFIC  VOLUME,  V 


TIME,  t 


JA-3 14543-20 A 


Figure  1 .  Pressure-specific  volume  paths  and  corresponding  particle  velocity 
histories  for  ideal  detonation  with  resolved  reaction  zone. 


The  SRI  large-scale  Lagrange  particle  velocity  gage  technique  used 
here  has  been  described  previously.^ The  basic  experimental  con¬ 
figuration  Is  shown  In  Figure  2.  The  HE  target  Is  a  right  cylinder 
typically  0.2-0. 3  m  In  diameter  by  about  0.13  ra  In  height.  A  plane-wave 
lens  and  appropriate  HE  driver  are  used  to  Initiate  steady  planar 
detonation  In  the  target  HE.  A  number  of  electromagnetic  particle 
velocity  (EPV)  gages, described  below,  are  embedded  at  various  depths 


(Lagrange  coordinates)  in  the  target  to  record  the  particle  velocities 
at  these  Lagrange  coordinates.  Expendable  Helmholtz  colls  wound  on 
plywood  frames,  and  powered  by  a  capacitive  discharge  system  with  a  time 
constant  of  several  milliseconds,  provide  a  uniform  uniaxial  magnetic 
Induction  of  400-500  gauss  In  the  region  of  the  EPV  gages.  The  magnetic 
induction  Is  essentially  constant  during  the  measurement  Interval,  about 
10  ps.  The  configuration  In  Figure  2  is  termed  large-scale  In  com¬ 
parison  to  Lagrange  gage  gas  gun  experiments  performed  In  other  labora¬ 
tories,  typically  with  0.05-m-diaraeter  targets,  to  study  shock,  initia¬ 
tion  of  detonation  for  shorter  times. 


j  HELMHOLTZ  COILS  ON  PLYWOOD  FRAME  £ 

MA-330I-1C 

Figure  2,  Configuration  of  SRI  large-scale  Lagrange 
particle  velocity  gage  experiments. 

Heavy  clashes  in  the  HE  target  indicate  the 
active  elements  of  the  embedded  particle 
velocity  gages. 

The  principle  of  the  EPV  gage  is  that  an  electrical  conductor 
moving  In  a  fixed  magnetic  field  will  generate  a  motional  electromotive 
force  (EMF)  proportional  to  Its  velocity.  If  the  conductor,  motion,  and 


19 


magnetic  induction  are  each  uniaxial  and  mutually  orthogonal,  as  in  the 
present  case,  then  by  Faraday's  law  of  Induction  for  moving  circuits 

E(t)  -  BJlu'(t)  (1) 

where  E  is  the  EMF,  t  is  time,  B  is  the  magnitude  of  the  magnetic  induc¬ 
tion,  X  is  the  length  of  the  conductor,  and  u'  is  its  velocity. 

The  EPV  gages  used  in  this  work  are  formed  from  aluminum  strips 
nominally  0.15  ram  thick  by  3  ram  wide.  The  strips  are  bent  into  rec¬ 
tilinear  U-shapes,  and  ten  are  mounted  on  a  standard  linen  phenolic  gage 
block  for  embedding  in  a  target  as  shown  in  Figure  3.  The  crossbar  of 
the  U  is  the  active  element  of  the  gage  and  is  nominally  25.4  mm  long. 
The  sides  of  the  U,  electrical  leads  that  carry  the  signal  out  the  back 
of  the  target,  do  not  contribute  to  the  signal  because  they  are  strictly 
parallel  to  the  direction  of  particle  motion  in  the  target.  The  gages 
are  anodized  and  sealed  to  provide  an  electrical  insulating  layer  5- 
10  pm  thick.  Because  the  active  elements  of  the  EPV  gages  are  only  0.15 
mm  thick  in  the  wave  propagation  direction  and  because  of  the  reasonable 
Impedance  match  between  the  aluminum  gage  material  and  the  HE  detonation 
products,  the  gages  reach  close  mechanical  equilibrium  with  instan¬ 
taneous  changes  in  the  surrounding  flow  in  the  products  in  less  than 
0.1  ps.  Therefore  u'(t)  is  taken  equal  to  the  local  particle  in  the 
reactive  flow,  u(t),  everywhere  except  at  the  spike  state,  which 
generally  cannot  be  fully  resolved  with  the  present  EPV  gages. 

A  high  frequency  recording  circuit,  consisting  of  RG  213  or  RG  58 
coaxial  cables  and  a  51-Q  vlewlng-and-termlnatlon  resistor  in  series 
with  the  EPV  gage,  monitors  each  gage.  Tektronix  7000  series  oscillo¬ 
scopes,  time  correlated  by  a  repetitive  beam-blanking  pulse  simul¬ 
taneously  applied  to  each  oscilloscope,  record  the  potential  drops 
across  the  viewing  resistors.  Typical  peak  potential  drops  are  2  V. 

E(t)  is  determined  from  the  potential  drop  across  the  viewing  resistor, 
and  u'(t).  Interpreted  as  u(t),  is  then  determined  from  Equation  (1). 


21 


Our  standard  ten-gage  gage  block  (Figure  3)  was  used  in  this 
work.  The  gage  numbering  system  and  dimensions  In  the  wave  propagation 
direction,  that  Is,  along  the  Lagrange  coordinate  axis,  are  shown  in 
Figure  4.  The  gages  are  numbered  sequentially  from  one  end  of  the  block 
to  the  other.  The  ten  gages  consist  of  two  symmetrically  placed  5-gage 
sets,  1  through  5  and  6  through  10.  Thus  there  are  two  replicate  gages, 
one  from  each  5-gage  set  In  each  of  five  planes  (Lagrange  positions) 
orthogonal  to  the  detonation  propagation  direction.  Adjacent  gages  are 
separated  by  12.7  ±  0.1  mm  In  the  propagation  direction  and  by  12.0  i 
0.2  mm,  center-to-center,  in  the  lateral  direction. 

2.2  GUINSY  UNIAXIAL  FLOW  LAGRANGE  ANALYSIS  CALCULATIONS 

The  Lagrange  analysis  for  uniaxial  flow  is  a  method  for  determining 
p,  V,  and  e  (specific  Internal  energy)  histories  and  cross  plots  from 
the  experimentally  measured  u(t)  data  obtained  at  different  Lagrange 
coordinates  h.  The  analysis  numerically  integrates  the  differential 
form  of  the  conservation  relations  for  momentum,  mass,  and  energy  In 
uniaxial  one-dlmenslonal  (ID)  flow.  These  relations.  In  a  form  ready 
for  integration  from  particle  velocity  gage  data,  may  be  written  as 


where  t  »  Lagrange  time  (the  time  a  particle  h  enters  the  flow),  H(t)  Is 
the  shock  path  In  the  h-t  plane,  and  the  subscript  0  denotes  the  undis¬ 
turbed  state  in  the  unreacted  explosive  ahead  of  the  shock. 

The  SRI  GUINSY  code^  solves  these  relations  using  the  experimental 
u(t)  data  at  various  h  to  approximate  the  required  spatial  and  time 
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derivatives  of  particle  velocity.  GUINSY  results  based  on  particle 
velocity  histories  for  steady  detonation  have  been  shown  to  be  accurate 
by  comparison  with  analytical  solutions.^ 

Because  the  uniaxial  flow  Lagrange  analysis  is  based  on  solution  of 
the  fundamental  conservation  relations  without  recourse  to  constitutive 
relations,  it  is  general  and  applies,  without  restrictions,  to  any 
reactive  or  nonreactive  uniaxial  strain  flow.  Because  the  measurements 
used  to  obtain  the  solutions  are  made  within  the  unperturbed  flow  rather 
than  at  a  free  surface  or  in  a  witness  material,  no  unfolding  procedures 
are  necessary  to  relate  the  results  to  the  basic  detonation  process. 

The  Lagrange  gage  and  analysis  technique  is  the  most  direct  and  accurate 
method  we  know  to  determine  CJ  parameters. 

2.3  TIGER  THERMODYNAMIC  EQUILIBRIUM  CALCULATIONS 

O 

TIGER  is  a  thermodynamic  equilibrium  code  developed  at  SRI  for 
calculating  the  equilibrium  thermodynamic  properties  of  a  heterogeneous 
system  of  known  atomic  composition  governed  by  an  arbitrary  equation  of 
state.  Assuming  thermodynamic  equilibrium  (mechanical,  thermal,  and 
chemical)  for  the  product  species  and  conservation  of  numbers  of  atoms 
of  each  element  between  reactants  and  products,  TIGER  calculates  the 
product  composition  (specific  mole  numbers  of  each  allowed  species)  and 
the  remaining  thermodynamic  state  variables  if  two  state  variables  such 
as  p  and  v  are  specified.  Furthermore,  under  the  same  assumptions, 

TIGER  calculates  the  thermodynamic  state  variables  and  the  product 
composition  along  paths  on  the  equation  of  state  surface  corresponding 
to  well  defined  thermodynamic  processes  such  as  shocks,  CJ  detonation, 
and  Isentroplc  expansion.  Thus,  in  this  work  we  are  able  to  estimate 
Hugonlots,  CJ  parameters,  and  Isentropes  through  specified  p-v  states 
using  TIGER. 

Two  equations  of  state  for  gaseous  species  are  commonly  used  in 
TIGER  caculatlons:  BKW  (Becker-Kistlakowsky-Wilson)  and  JCZ3  (Jacobs- 
Cowperthwalte-Zwlslei ) .  These  are  both  nonideal  gaseous  p-v-T-n^^  equa¬ 
tions  of  state,  where  T  is  temperature  and  the  nj^  are  the  mole  numbers 
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of  the  various  species  present  In  the  products.  Such  EOSs  that  contain 
the  product  composition  as  Independent  variables  are  often  said  to  have 
explicit  chemistry.  Both  BKW  and  JCZ3  represent  the  state  of  the  art 
In  such  EOSs.  JCZ3,  however,  has  been  formulated  to  predict  a  lattice 
pressure  at  0  K  and  BKW  has  not.  In  practice,  there  are  quantitative 
differences,  often  small,  between  their  predicted  values  for  detonation 
parameters.  We  find  that  TIGER/JCZ3  predictions  are  usually  closer  to 
experimental  values,  especially  for  detonation  velocity.  Further 
examples  and  discussion  are  given  by  Finger  et  al.^® 

In  this  work  TIGER  calculations  using  both  BKW  and  JCZ3  are  used  to 
model  and  extend  the  results  obtained  for  each  of  the  two  study  mater¬ 
ials  with  the  Lagrange  gage  and  analysis  techniques.  Procedures  vary 
considerably  In  the  two  cases  and  are  described  In  Sections  3  and  4. 


2.4  JWL  CJ  ISENTROPE  FITTING  PROCEDURE 

JWL  Is  an  empirical  detonation  products  EOS  without  explicit 
chemistry.  The  form  of  the  EOS  Is 

-R,V  -R-,V  „ 

A/1  \  1  D/1  “  \  2  ,  0)  E 

p  -  A(1  -  )e  +  B(1  -  )e  +  ^  (5) 

and  the  corresponding  JWL  form  for  an  Isentrope  Is 

-(w+l) 


-R, V  -R-V 
Pg=  Ae  +  Be  +  CV 


(6) 


where  E  Is  specific  energy /volume  (E  *  s/vg)  measured  with  respect  to 
the  energy  level  of  fully  expanded  (v  =  Infinity)  products;  V  Is  rela¬ 
tive  volume  (V  =  v/vq)  and  Is  dimensionless;  A,  B,  and  C  are  the  linear 
JWL  coefficients,  with  C  being  an  Integration  constant  assigned  the 
appropriate  value  to  specify  the  CJ  Isentrope;  and  R^ ,  R2,  and  w  are  the 
nonlinear  JWL  coefficients. 

JWL  Is  an  E(p,V)  EOS  of  the  Mle-Grunelsen  type  with  the  reference 
curve  being  the  CJ  Isentrope,  Equation  (6),  and  the  Grunelsen  parameter 
being  the  constant  u).  The  form  of  the  JWL  CJ  Isentrope  Is  such  that. 


for  large  V  (In  practice,  V  >  6-8),  the  exponential  terms  become 
completely  negligible.  Therefore,  for  large  V  the  JWL  CJ  Isentrope 
takes  the  form  of  the  polytropic  gas  CJ  Isentrope,  pV^  =  const,  where 
the  polytropic  constant  y  Is  the  ratio  of  specific  heats  at  constant 
pressure  and  volume  and  from  Equation  (6),  y  “(“■*■  !)• 

Thus,  the  value  of  the  nonlinear  coefficient  u)  strongly  affects  the 
detonation  product  behavior  predicted  by  JWL  both  for  large  V  on  the  CJ 
Isentrope  (the  low  pressure  polytropic  limit)  and  for  all  V  off  the  CJ 
Isentrope  (the  constant  Grunelsen  parameter).  Along  the  CJ  Isentrope, 
the  Rj^  term  often  dominates  behavior  for  V^j  <  V  <  1,  and  all  terms  are 
significant  for  1  <  V  <  ~5. 

The  present  form  of  JWL  was  developed  by  Lee  et.  al.^^  to 
accurately  predict  the  metal  pushing  behavior  of  HEs  In  spherical  and 
cylindrical  geometries  to  large  V.  The  six  parameters  are  usually  fit 
to  reproduce  the  results  of  cylinder  tests^^*^^  to  expansions  approach¬ 
ing  10  as  well  as  to  satisfy  certain  thermodynamic,  CJ  state,  and  avail¬ 
able  energy  constraints  described  In  reference  11.  General  properties 

of  JWL  are  reviewed  by  Flckett  and  Davis, ^  and  JWL  parameters  for  a 
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number  of  HEs  are  given  by  Dobratz. 


The  objective  of  our  work  was  to  develop  Initial  JWL  parameters  for 
certain  new  HEs  so  that  hydrocode  calculations  could  be  performed  to 
evaluate  the  behavior  of  these  HEs  In  various  NWE  simulators.  In  any 
one  application,  the  detonation  products  may  undergo  different  amounts 
of  volume  expansion  and  experience  different  energy  exchange  mechanisms 
with  the  surroundings  than  they  do  In  another  or  In  the  cylinder  test. 
Therefore,  the  present  characterization  results  were  not  based  on  the 
results  of  cylinder  tests  or  other  experimental  or  computational  simula¬ 
tions  of  a  specific  potential  DNA  application.  Rather,  to  generate 
broadly  applicable  JWL  parameters,  we  performed  experiments  and  computa¬ 
tions  to  obtain  an  estimate  of  the  general  CJ  isentrope  over  as  large  a 


range  of  volume  expansions  as  possible  and  then  developed  a  procedure, 
described  below,  to  uniquely  fit  the  JWL  form  to  this  whole  estimated  CJ 
"basis  Isentrope." 
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If  more  accurate  JWL  parameters  are  desired  for  a  specific  applica¬ 
tion,  the  CJ  basis  Isentrope  estimate  can  be  Improved  In  the  range  of 
Interest,  for  example,  by  performing  a  relevant  simulation  experiment  to 
constrain  the  specification  of  Isentroplc  release  states  in  this 
range.  The  fitting  techniques  developed  here  can  then  be  readily 
adapted  to  the  new  characterization  range  to  uniquely  determine  JWL  EOS 
parameters  optimized  to  the  application. 

Fitting  Procedure 

The  objective  of  the  fitting  procedure  Is  to  determine  values  for 
the  six  JWL  EOS  coefficients:  A,  B,  C,  R^,  R2»  and  o)  and  to  specify  the 
four  CJ  state  coefficients  with  which  they  are  mathematically  consis¬ 
tent:  p,  0,  r,  and  Eq. 

r  Is  the  adiabatic  exponent  at  the  CJ  state  defined  as  the 
logarithmic  derivative  (d  In  pg/d  In  V)cj.  It  serves  to  specify  Vqj, 
the  relative  volume  of  the  products  at  the  CJ  state  since,  as  a  result 
of  the  tangency  of  the  CJ  Isentrope  and  the  products  Hugoniot  at  the  CJ 
state,  =  r/(r  +  1). 

Eq  is  defined  by  the  equation  Eqj  =  Eg  +  (pgj/2)(l  -  Vqj).  In  the 
cylinder  test  approach.  Eg,  which  Is  closely  related  to  the  maximum 
useful  work-potential  of  an  HE  in  a  detonation,^  is  chosen  consistent 
with  available  chemical  energy  as  determined  either  from  detonation 
calorimetry  or  from  thermodynamic  calculations.  Through  Its  defining 
equation.  Eg  then  specifies  E^j  and  provides  one  of  the  constraints  used 
to  determine  the  JWL  parameters. In  the  present  approach.  Eg  is  also 
consistent  with  thermodynamic  calculations,  but  is  a  derived  quantity 
determined  from  the  JWL  parameters  as  described  below. 

The  basis  for  the  present  fitting  procedure  Is  the  set  of  points  In 
the  p-V  plane  defining  the  estimated  CJ  Isentrope  between  '^CJ  and  some 
maximum  relative  volume,  Vj.  Call  this  "basis"  Isentrope  defined 

for  Vgj  <  V  <  Vf.  (Values  of  a  variable  pertaining  to  the  basis  Isen¬ 
trope  are  now  written  without  superscript;  values  of  that  variable 
pertaining  to  the  JWL  fit  are  written  with  the  superscript  J.  The 


subscript  S  denotes  the  CJ  Isentrope.)  Let  A  =  Eg(V^j)  -  E3(Vj)  and 
note  that  we  can  determine  A  from  the  basis  Isentrope  using  the  first 
law  of  thermodynamics  for  an  adiabatic  process 


PgCV)  dV 
f 


In  fitting  the  JWL  CJ  Isentrope  form.  Equation  (6),  to  the  basis  CJ 
Isentrope,  we  Impose  three  strict  mathematical  constraints.  First, 


p/(v..) 


CJ^ 


(7) 


The  CJ  pressure-relative  volume  points  on  the  basis  and  JWL  CJ  Isen- 
tropes  coincide.  Second, 


Ps^^CJ^ 


Ps  ('^CJ> 


1  -  V, 


CJ 


d  V 


(8) 


The  JWL  CJ  Isentrope  and  the  basis  Rayleigh  line  (locus  of  material 
states  in  shock  and  steady-state  reaction  zone)  are,  by  Equations  (7) 
and  (8),  tangent  in  the  p-V  plane  at  the  CJ  state.  Third, 


(9) 


The  volumetric  Internal  energy  differences  on  the  JWL  and  the  basis  CJ 
Isentropes  between  the  CJ  relative  volume  and  the  final  "basis"  relative 
volume  are  equal - 

In  addition  to  these  mathematical  constraints,  following  Lee  et 
al.,^^  we  Impose  two  further  restrictions  on  allowed  fits.  First,  w 
must  be  In  the  range  0.2  <  w  <  0.4,  and  second,  the  linear  coefficients 
must  be  positive. 

The  following  solution  procedure  provides  unique  values  for  the  six 
JWL  EOS  parameters;  these  values  constitute  a  fit  that  closely  repro¬ 
duces  the  basis  CJ  Isentrope  and  precisely  satisfies  the  three 
constraints  above  subject  to  the  two  cited  restrictions.  First,  we 


select  a  threshold  relative  volume  V^.  in  the  range  10  <  such 

that  for  V  >  V^,  pgCV)  Is  well  represented  as  a  constant  y  gas  with 

1.2  <  y  <  1.4.  Since  >  10,  we  know  that  for  V  >  V  the  exponential 

terras  in  Pg'^(V),  Equation  (6),  are  completely  negligible.  Therefore, 

for  V  <  V  <  V  ,  we  perform  a  least  squares  fit  of  the  final  term  for 
I  ^  ^ 

Pg  (V)  in  Equation  (6)  to  the  basis  Isentrope  and  thereby  determine  C 
and  Cl).  The  resulting  value  for  co  satisfies  the  first  restriction  stated 
above,  assuring  a  reasonable  value  for  the  low  pressure  adiabatic 
exponent. 

To  determine  the  remaining  four  JWL  EOS  parameters,  we  next  use 
constraining  Equations  (7)  and  (8)  to  express  the  two  undetermined 
linear  coefficients,  A  and  B,  as  functions  of  the  two  unknown  nonlinear 
coefficients  Ri  and  R2»  and  known  quantities.  We  use  these  expressions 
and  the  third  constraint.  Equation  (9),  to  obtain  an  implicit  equation 
for  R]^  as  a  function  of  R2*  this  step  we  make  use  of  the  fact  that 
E3'^(V£)  is  a  calculable  number  because  and  C  and  w  have  been 

evaluated.  To  obtain  an  additional  relation  between  R^  and  R2  com¬ 
pleting  the  system  of  equations  and  allowing  us  to  determine  the  four 
unknown  coefficients,  we  require  that  the  RMS  deviation  between  the 
basis  and  JWL  CJ  Isentropes  be  a  minimum. 

The  solution  procedure  described  above  determines  the  six  JWL  EOS 
coefficients.  The  only  step  that  is  not  objective  is  the  selection  of 
V^.  The  fitting  procedure  is  completed  by  calculating  Eq  from  its 
defining  equation  and  verifying  that  all  the  linear  coefficients  are 
positive.  Negative  coefficients  were  not  encountered  wl'h  this  fitting 
procedure  and  the  basis  CJ  Isentropes  used  in  this  work,  but  they  can 
arise.  If  they  do,  appropriate  modifications  must  be  made  either  to  the 
basis  isentrope  or  to  the  fitting  procedure. 

Discussion 

Although  the  present  fitting  procedure  provides  unique  values  for 
the  JWL  EOS  coefficients  and  closely  reproduces  the  basis  CJ  isentrope, 
other  equally  valid  procedures  producing  different  unique  fits  are 


possible.  For  example,  another  measure  of  goodness  of  fit  between 
pressures  on  the  basis  and  JWL  CJ  Isentropes  could  be  Invoked  or,  more 
fundamentally,  a  different  property  of  the  basis  CJ  Isentrope  such  as 
the  adiabatic  exponent  could  be  matched.  In  the  present  approach  the 
Y  -  V  behavior  Is  not  controlled;  It  Is  a  result  of  the  fitting  proce¬ 
dure.  The  motivation  for  the  present  approach  Is  our  belief  that  the 
locus  of  the  CJ  Isentrope  Is  the  most  important  parameter  for  estimating 
the  performance  of  an  HE  In  a  simulator. 

The  first  two  constraints.  Equations  (7)  and  (8),  are  the  same  as 
those  used  in  fitting  the  JWL  EOS  parameters  from  cylinder  test 
data.^^  The  energy  constraint.  Equation  (9),  however.  Is  different. 

The  resulting  Eq  value  In  our  procedure  is  completely  consistent  with 
the  JWL  EOS  form  and.  In  the  cases  to  which  It  Is  applied  here,  with 
thermodynamic  equilibrium  calculations  to  V^.  However,  we  do  not  Invoke 
the  condition  that  Eg  equal  to  the  available  chemical  energy 
estimated  from  the  product  composition  at  some  point  in  the  equilibrium 
calculations  and  from  standard  heats  of  formation  of  the  products  and 
reactants.  In  fact,  our  Eq  tnay  be  less  or  more  than  the  available 
chemical  energy.  This  difference  between  the  present  technique  and  the 
cylinder  test  techniques  for  fitting  JWL  parameters  reflects  the 
behavior  of  the  products  between  =  200-1000  and  V  »  »,  occurlng 
mostly  at  p  <  1  atm,  and  Is  judged  to  be  of  no  concern  In  DNA  NWE  simu¬ 
lator  hydrocode  calculations  with  JWL  fits  from  either  method. 

The  fitting  procedure  developed  here  Is  believed  to  be  general  and 
useful  In  a  number  of  situations.  The  optimization  approach  generates 
unique  values  for  all  the  JWL  EOS  parameters  and  can  be  adapted  to  other 
characterization  data  bases  such  as  cylinder  test  results.  The  specific 
method  used  here  with  the  CJ  Isentrope  is  Independent  of  the  source  of 
the  basis  CJ  Isentrope  and  is  rapid  and  straightforward  to  apply.  Thus, 
It  provides  a  convenient  means  for  transforming  the  results  of  thermo¬ 
dynamic  calculations  based  on  nonideal  EOSs  with  explicit  chemistry  into 
the  computationally  preferable  nonexpllclt  JWL  form.  The  resulting  JWL 
EOS  Is  as  consistent  as  possible  with  the  explicit  EOS  In  the  Important 


region  on  and  near  the  CJ  Isentrope  and  provides  a  good  Initial  computa¬ 
tionally  convenient  characterization  for  a  new  HE  In  the  absence  of 
experimental  data. 


SECTION  3 


NA/NP  CHARACTERIZATION 

Initial  JWL  coefficients  were  developed  for  a  representative  member 
of  the  NA-NP  family,  liquid  explosives  similar  In  performance  to  nltro- 
methane  (NM).  The  specific  composition  selected  for  study  was  a  blend 
Identified  In  prior  commercial  studles^^”^^  as  offering  a  favorable  com¬ 
bination  of  Initiation,  performance,  cost,  and  handling  properties. 
Throughout  this  report,  we  use  the  term  NA/NP  to  refer  to  this  baseline 
blend,  the  present  study  material. 

The  tasks  performed  In  developing  JWL  coefficients  for  NA/NP 
Included  development  of  Initial  properties  and  handling  procedures  for 
the  study  material;  Lagrange  gage  detonation  and  shock  Initiation 
experiments  and  Lagrange  analysis  to  determine  peak,  CJ,  and  Initial  CJ 
release  Isentrope  states;  TIGER  calculations  to  extend  the  CJ  release 
Isentrope  to  large  volumes;  and  determination  of  JWL  coefficients.  The 
following  subsections  present  the  results  of  these  tasks  and  conclude 
with  a  discussion. 

3.1  DESCRIPTION  OF  MATERIAL 
Background 

NA/NP  belongs  to  a  family  of  liquid  solutions  consisting  of  nitric 
acid,  nltropropane  (NP),  and  water.  It  Is  convenient  for  this  report, 
and  In  practice,  to  treat  the  nitric  acid  and  the  water  as  a  single  com¬ 
ponent,  an  aqueous  solution  of  nitric  acid  referred  to  here  as  NA.  The 
properties  of  the  NA-NP  family  relevant  to  Its  potential  use  as  an  HE 
were  Initially  examined  by  Mallory  under  contract  to  Troclno  and 
Associates. Troclno  later  compiled  Mallory's  results  as  well  as 
practical  data  and  handling  requirements  for  NA  and  NP  In  a  very  useful 
reference  on  NA-NP  as  an  HE  for  military  and  commercial  applications.^^ 
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Mallory’s  investigations,  briefly  summarized  below,  are  the  only 
research  on  NA-NP  of  which  we  are  aware  other  than  this  work. 


For  NA-NP  blends  o^qrgen-balanced  to  CO  or  CO2,  Mallory  first 

examined  the  effect  of  NA  concentration,  over  the  practical  range  from 

60  to  90  wt%,  on  properties  including  gap  test  sensitivity,  plate  dent, 
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detonation  velocity,  and  NP  separation  temperature  (solubility  limit). 
This  work  indicated  that  the  CO2  balanced  compositions  with  acid  concen¬ 
trations  near  or  slightly  greater  than  65%  were  stable  solutions  at  nor¬ 
mal  temperatures  and  had  detonation  velocities  and  gap  sensitivities 
slightly  greater  than  neat  NM.  Increasing  acid  concentrations  led  to 
lower  separation  temperatures,  greater  detonation  velocities,  and  much 
greater  gap  sensitivities,  with  the  last  apparently  exceeding  that  of 
tetryl  for  some  blends  with  NA  concentrations  of  90%.  A  major  conclu¬ 
sion  was  that  the  abundant,  inexpensive,  manageable  NAs  with  relatively 
low  acid  concentrations  near  65%,  in  CO2  balanced  mixtures,  produced  NA- 
NPs  with  a  collection  of  good  HE  properties. 

In  his  second  report, Mallory  focused  on  the  low  concentration 
NAs  and  examined  NA-NP  properties  including  cap  and  booster  sensitivity, 
mixture  variation  effects,  calculated  CJ  pressure  and  detonation  velo¬ 
city  based  on  the  Kamlet  and  Jacobs  Ruby  code  approximation  method, 
and  performance  in  shaped  charges.  This  work  identified  Baume  42  NA 
(concentration  >•  67.18  wt%)  and  NP,  CO2  balanced  by  blending  2.24  parts 
NA  and  1  part  NP  by  volume,  as  the  most  promising  from  all  aspects. 

This  blend,  described  below,  is  NA/NP,  the  present  study  material. 

The  results  of  the  detonation  experiments  and  calculations  per¬ 
formed  here  are  in  good  agreement  with  Mallory's  results.  Specific 
comparisons  are  given  as  the  results  are  presented. 

Present  Study  Material 

The  present  study  material  is  Baume  42  NA  and  NP  blended  to  be  oxy¬ 
gen  balanced  to  CO2.  Note  that  our  specification  of  the  composition  of 
this  blend  is  slightly  different  from  Mallory's.  Below  we  discuss 


Baume  42  nitric  acid,  the  NP  uaed  In  this  work,  the  NA/NP  oxygen  balance 
condition  and  resulting  composition,  the  mixture  precision  for  our 
NA/NP,  the  Initial  density  of  our  NA/NP,  and  our  calculated  heats  of 
formation  and  detonation  for  NA/NP. 

Baume  42  nitric  acid  Is  an  Industrial  grade  aqueous  solution  of 
nitric  acid  and  Is  available  from  several  manufacturers.  In  a  Du  Pont 
technical  bulletin, It  Is  termed  weak  nitric  acid,  technical;  at 
15.6“C  (60“F)  it  Is  stated  to  be  (nominally)  67.18  wt%  nitric  acid  with 
a  specific  gravity  of  1.4087  and  a  temperature  correction  factor  of 
0.00084/®F.  (Because  the  Baume  scale  Is  a  measure  of  density,  the 
correspondence  between  an  actual  Baume  reading  and  the  associated  acid 
concentration  Is  a  function  of  temperature;  however.  In  this  case  the 
Baume  42  value  Is  being  used  simply  as  a  label  for  NA  of  concentration 
67.18  wt%,  and  this  association  la  therefore  valid  at  all  temperatures.) 

On  the  basis  that  Baume  42  NA  Is  0.6718  by  weight  pure  nitric  acid 
and  0.3282  water,  we  calculated  the  following  molecular  composition  and 
heat  of  formation.  Using  the  1959  atomic  weight  scale  (0  =  16)  In 
conformity  with  Mallory's  procedures,  there  are  1.7088  molecules  of 
water  per  molecule  of  pure  nitric  acid  In  Baume  42  NA.  The  heat  of 
formation.  Including  heat  of  dilution.  Is  -162.63  kcal/raole.  In  agree¬ 
ment  with  Mallory's  value  less  the  heat  of  dilution. 

Baume  42  NA  appears  to  be  one  of  the  least  expensive  and  most 
easily  handled  forms  of  nitric  acid.  We  encountered  no  difficulties  In 
storing  or  handling  It,  although  NAs  of  higher  or  lower  concentration 
might  present  additional  problems  related  to  toxicity  or  activity 
respectively.  In  the  Lagrange  gage  experiments  at  our  remote  site,  it 
was  convenient  to  perform  the  NA  handling  operations  outdoors  while 
wearing  conventional  protective  gear  Including  gloves,  face  shields,  and 
disposable  overgarments.  Over  the  preparation  periods  for  our  experi¬ 
ments  (24  hours  or  less),  we  detected  no  reaction  between  the  NA  and  our 
anodized  aluminum  gages  or  the  polyethylene  vessels  we  used  for  tem¬ 
porary  NA  storage  or  weighing  containers. 


The  NP  used  here  was  furnished  by  the  International  Minerals  and 
Chemical  Corporation,  IMC  Chemical  Group,  666  Garland  Place,  Des 
Plaines,  IL  60016  (sample  number  DP-8562-TD).  It  Is  described  as  a 
mixture  of  the  two  stable  Isomers  for  each  of  which  properties  are  given 
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In  two  technical  notes.  ♦  We  attributed  the  following  properties  to 
the  mixture  by  averaging  values  for  the  Isomers.  Specific  gravity  at 
ZS^C  =  0.992.  Density  correctlon/“C  =  0.0011  Mg/m^“C.  Heat  of  forma¬ 
tion  =  -41.7  kcal/mole.  (In  some  TIGER  calculations,  we  used  the  value 
for  the  Isomer  NP-2,  namely,  -43.2  kcal/mole,  estimated  to  have  a  negli¬ 
gible  effect  on  the  results.) 

We  found  NP  to  be  easy  to  handle.  In  our  preshot  preparations 
performed  outdoors,  no  special  procedures  were  used  other  than  prevent¬ 
ing  contact  between  the  NP  and  those  polymers  for  which  It  Is  a  strong 
solvent.  We  stored  It  In  polyethylene  containers  and  used  PVC  or  poly¬ 
ethylene  foil  or  wax  In  our  experimental  assemblies  to  Isolate  It  from 
components  It  could  dissolve. 

NA/NP  Is  said  to  be  oxygen  balanced  to  CO2  In  that  the  reactant 
composition  Is  chosen  to  produce  only  the  stable  detonation  products 
CO2,  H2O,  and  N2.  For  NA/NP,  including  the  water  In  the  NA,  the  CO2 
oxygen  balance  equation  Is 

3HNO3  +  3(1.7088)H20  +  C3H7NO2  3C02  +  [5  +  3(1 .7088) ]H20  +  2N2 

Given  the  CO2  oxygen  balance  equation,  the  molecular  weights  of 
pure  nitric  acid  and  NP,  and  the  concentration  of  pure  nitric  acid  In 
Baume  42  NA,  we  can  calculate  the  composition  for  NA/NP.  The  result  Is 
that  NA/NP  Is  3.158  parts  NA  to  1  part  NP,  by  weight.  Alternatively, 
the  weight  fraction  of  NA  In  NA/NP  Is  0.7595. 

We  prefer  to  specify  the  NA/NP  composition  by  weight  because  these 
values  are  Independent  of  temperature.  The  composition  by  volume,  as 
given  by  Mallory,  has  a  slight  temperature  dependence.  For  example,  at 
the  temperature  of  our  experiments,  29®C,  the  composition  by  volume, 
calculated  from  the  density  and  specific  gravity  Information  cited  above 
for  NA  and  NP,  Is  2.251  parts  NA  to  1  part  NP. 
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We  mixed  our  NA/NP  from  preweighed  stocks  of  NA  and  NP  having  the 
weight  ratios  given  above*  The  stocks  were  weighed  on  an  electronic 
balance  with  a  precision  of  0.1  g.  The  balance  was  calibrated  after 
each  weighing  with  a  collection  of  standard  weights  chosen  to  have  an 
aggregate  weight  about  equal  to  the  stock  plus  container  weight.  These 
aggregate  weights  were  between  1  and  5  kg.  We  estimate  that  the  weights 
of  each  component  stock  were  controlled  to  accuracies  better  than  0.1% 
by  this  method.  Therefore,  the  present  NA/NP  study  material  Is  judged 
to  be  well  within  0.5%  of  the  calculated  oxygen  balance  condition.  (The 
densities  of  the  NA  and  the  NP  were  checked  and  agreed  with  the  nominal 
values  given  previously  within  0.3%.) 

We  measured  the  density  of  the  mixed  NA/NP  In  the  Lagrange  gage 
experiments  to  be  1.252  ±  0.005  Mg/m^  at  29'’C.  We  used  this  value  In 
our  calculations  although  the  last  significant  figure  is  not  reliable. 
From  the  nominal  densities,  specific  gravities,  and  correction  factors 
for  NA  and  NP,  and  Ignoring  mixing  effects  on  volume,  we  calculate  a 
density  of  1.260  Mg/m^  at  this  temperature,  in  good  agreement  with  our 
measurement*  If  we  use  the  measured  rather  than  nominal  densities  for 
NA  and  NP  and  apply  the  temperature  corrections,  the  calculated  density 
is  1.256  Mg/ra'^  in  even  better  agreement  with  the  measurement.  From  this 
we  conclude  that  the  present  solution  of  Baume  42  NA  and  NP  is  essen¬ 
tially  ideal,  although  Mallory  reports  a  detectable  volume  interaction 
on  mixing.  The  calculated  NA/NP  density  given  by  Mallory  is  1.272 
Mg/m^,  which  is  higher  than  our  calculations  but  pertains  to  NA/NP  at  a 
lower  temperature. 

From  the  heats  of  formation  for  reactants  given  previously,  we 
calculate  a  heat  of  formation  of  -176.53  kcal/g  and  a  maximum  heat  of 
detonation  of  911  kcal/g.  Because  NA/NP  is  oxygen  balanced  to  CO2,  the 
detonation  product  composition  is  nearly  constant  throughout  the  product 
expansion;  thus  the  maximum  heat  of  detonation  is  close  to  the  effective 
heat  of  detonation  of  NA/NP. 


3.2  LAGRANGE  GAGE  EXPERIMENTS 


Three  large-scale  multiple  Lagrange  particle  velocity  gage  experi¬ 
ments  were  performed  on  NA/NP.  Two  were  detonation  experiments,  and  the 
third  was  a  shock  Initiation  experiment.  The  targets,  Initiation  sys¬ 
tems,  and  results  for  these  experiments  are  described  below. 

NA/NP  Target  Chambers 

The  target  chamber  design  used  In  all  the  NA/NP  experiments  Is 
shown  schematically  In  Figure  5.  The  chambers  were  fabricated  from 
Pyrex  beakers  and  were  approximately  0.19  m  In  Inside  diameter  by  0.12  m 
high.  Precision  ground  parallel  surfaces  on  the  tops  and  bottoms 
ensured  that  the  gage  elements  were  parallel  to  the  HE  driver  face  to 
within  0.02  mm.  The  joints  and  the  tops  of  the  linen  phenolic  gage 
blocks  were  sealed  with  wax  to  protect  against  the  solvent  effects  of 
the  NP.  A  1 .42-mm-thlck  PMMA  plate,  protected  on  the  Inside  by  0.01-mm- 
thlck  PVC  foil,  separated  the  NA/NP  from  the  driver  HE.  Each  chamber 
contained  a  fill  tube  and  a  bleed  tube  to  remove  air  bubbles.  A  photo¬ 
graph  of  one  of  the  chambers  Is  shown  In  Figure  6. 

The  filling  procedure  was  to  mix  the  preweighed  stocks  of  NA  and  NP 
at  the  time  of  the  experiment.  The  target  chamber  was  filled,  and  the 
experiment  performed.  In  the  orientation  Indicated  In  Figure  5,  that  Is, 
with  the  gage  block  down.  Therefore,  It  was  critical  to  avoid  entrapped 
air.  Air  bubbles  were  removed  through  the  bleed  tube,  and  fluid  heads 
were  left  In  each  tube  to  account  for  possible  leaks,  though  none  were 
detected. 

The  EPV  gages  In  the  NA/NP  experiments  were  hard  anodized.  Anodi¬ 
zation  was  considered  particularly  Important  In  these  experiments 
because  NA/NP,  unlike  other  HEs  we  have  examined.  Is  an  electrical 
conductor  In  the  unreacted  state.  Because  of  difficulties  encountered 
by  the  vendor  In  making  good  connections  to  each  of  the  gages  premounted 
on  the  gage  block,  the  gages  had  to  be  anodized  and  stripped  twice 
before  satisfactory  coatings  were  obtained.  Consequently,  the  gages  In 
these  experiments  were  somewhat  thinner  and  more  embrittled  than 
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normal.  In  the  later  NFN  experiments,  a  special  holding  fixture  was 
built  and  the  gages  were  anodized  before  they  were  attached  to  the 
precision  gage  blocks.  No  difficulties  were  encountered  with  this 
system. 

Initiation  Systems 

The  NA/NP  experiments  were  performed  with  0.2-m-dlameter  (8  Inch) 

HE  drivers  and  P-80  plane  wave  lenses.  We  selected  the  driver  HE  as 
follows.  We  assumed  that,  for  rapid  and  reliable  Initiation  of  steady 
detonation  In  NA/NP,  the  pressure  amplitude  of  the  initial  shock  induced 
In  the  unreacted  NA/NP  by  the  driver  (ignoring  the  spike  in  the  driver) 
should  be  approximately  equal  to  the  predicted  CJ  pressure.  This  condi¬ 
tion  Is  discussed  further  In  regard  to  the  NPN  experiments  In  Section 
4.2. 

We  very  approximately  estimated  a  Hugonlot  for  unreacted  NA/NP  by 
scaling  the  Hugonlot  pressures  of  unreacted  NM  according  to  the  density 
difference  between  NA/NP  and  NM.  We  then  graphically  performed  an 
Impedance  mismatch  calculation  for  TNT  and  Baratol  drivers  using  their 
estimated  products  CJ  release  isentropes,  as  shown  In  Figure  7.  The 
shock  amplitude  estimated  for  a  TNT  driver  was  158  kbar.  (The  CJ  pres¬ 
sure  estimated  by  Mallory  with  the  Ruby  code  approximation  method  of 
Kamlet  and  Jacobs  was  152.84  kbar.^^)  We  therefore  selected  TNT  for  the 
driver  HE  for  the  detonation  experiments.  The  effect  of  the  1.42-mm- 
thlck  PMMA  plate  was  neglected  because  of  the  approximate  nature  of 
these  calculations. 

By  the  same  techniques,  a  Baratol  driver  was  estimated  to  Induce  a 
much  lower  shock  pressure,  well  below  100  kbar,  In  unreacted  NA/NP  (see 
Figure  7).  A  Baratol  driver  was  therefore  used  in  the  Initiation 
experiment  to  generate  Initial  information  on  the  shock-to-detonation 
transition  (SOT)  In  NA/NP. 

Experimental  Results 

Three  multiple  Lagrange  particle  velocity  gage  experiments  were 
performed  on  NA/NP.  Shots  1  and  2  were  detonation  experiments  to 
observe  the  shock  compression,  reaction  zone,  and  initial  adiabatic 
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Figure  7.  Pressure-particle  velocity  plots  for  estimating  initiation 
pressures  in  NA/NP  experiments. 

All  curves  estimated;  do  not  use  for  other  applications. 
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expansion  states  during  steady  planar  detonation  In  NA/NP.  Shot  3  was  a 
preliminary  SDT  experiment.  The  results  of  the  detonation  experiments 
are  considered  first. 


Both  detonation  experiments  provided  a  good  data  return.  However, 
because  Improved  Instrumentation  settings  were  used  In  the  second 
experiment,  only  data  from  that  experiment  were  selected  for  quantita¬ 
tive  analysis.  Data  were  obtained  from  all  gages  except  10,  which  was 
found  to  be  open  before  the  experiment.  The  oscillographs  from  the 
second  NA/NP  detonation  experiment  are  shown  in  Figure  8. 

Vertical  sensitivities  for  the  oscillographs  In  Figure  8  are  500 
mV/dlvlsion;  horizontal  sensitivities  are  1  ps/dlvlslon.  Short  gaps  In 
the  profiles  are  the  repetitive  beam  blanks  used  to  time-correlate  the 
records . 

For  the  first  several  microseconds  after  shock  arrival,  the  quality 
of  the  gage  signals  from  this  experiment  Is  the  best  that  we  have  seen 
with  this  technique  (see  discussion  in  Section  4.2).  The  records  con¬ 
tain  exceptionally  little  noise  either  before  or  after  shock  arrival. 

There  are  only  two  negative  features  In  these  data  from  a  measure¬ 
ment  point  of  view.  First,  some  of  the  records  (for  example,  those  from 
gages  2,  3,  and  8)  show  a  slight  ramping  before  shock  arrival.  The 
cause  Is  not  known,  but  may  have  to  do  with  the  conductivity  of  the 
unreacted  NA/NP.  Second,  exceptionally  large  perturbations  occur  on  a 
given  record  when  the  next  downstream  gage  enters  the  flow  about  2  ps 
later.  This  effect  Is  so  pronounced  that  it  determines  the  useful 
recording  time  of  the  EPV  gages  In  this  experiment.  The  cause  Is  again 
not  known,  but  may  be  partial  breakdown  of  the  anodizing. 

The  nine  particle  velocity  histories  calculated  from  the  oscillo¬ 
graphs  In  Figure  8  are  shown  in  a  composite  plot  In  Figure  9.  For 
clarity,  the  histories  from  each  half  of  the  target  are  also  shown 
separately  In  Figure  10.  The  structure  of  these  histories  Is  considered 
In  the  analysis  presented  In  the  next  subsection. 
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Figure  9.  Lagrange  particle  velocity  histories,  gages  1  -9,  second  NA/NP 
detonation  experiment. 
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Figure  10.  Lagrange  particle  velocity  histories  for  each  half  of  target, 
second  NA/NP  detonation  experiment. 
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The  SDT  experiment  produced  unusual  results.  The  particle  velocity 
histories  are  shovm  In  a  composite  plot  in  Figure  11  and  for  each  half 
of  the  target  In  Figure  12.  These  histories  clearly  are  not  repro¬ 
ducible  across  the  target,  perhaps  Indicating  nonuniform  or  marginal 
Initiation.  The  peak  particle  velocity  amplitudes  are  less  and  the 
structures  of  both  the  Individual  profiles  and  the  flow  evolution  are 
qualitatively  different  from  those  for  the  detonation  experiment.  The 
histories  from  gages  1-5  show  no  obvious  pattern;  those  from  gages  6-10 
seem  to  consist  of  a  fairly  steady  spike  followed  by  a  growing  pulse. 
Histories  6-10  could  be  argued  to  be  tending  toward  the  shape  of  the 
profiles  In  the  steady  detonation  experiments. 

The  SDT  experiment  Is  not  analyzed  further  here  because  the  data 
are  ambiguous  and  because  the  primary  objective  of  this  task  was  to 
generate  JWL  coefficients  for  NA/NP.  However,  the  results  Indicate  that 
the  SDT  process  In  NA/NP  may  have  unusual  features,  such  as  a  surpris¬ 
ingly  long  time  to  reach  steady  detonation  and  possibly  a  complex  hydro- 
dynamic  phenomenology.  These  features  could  be  significant  In  DNA  NWE 
simulator  applications  and  should  be  Investigated. 

3.3  LAGRANGE  ANALYSIS 

The  particle  velocity  histories  from  the  second  NA/NP  detonation 
experiment.  Figures  9  and  10,  were  compared  to  obtain  a  single  data  set 
for  Lagrange  analysis.  The  first  step  was  to  discard  the  history  from 
gage  1  because  this  record  Is  anomalously  higher  than  the  other  records, 
and  there  Is  no  redundant  gage  history  to  corroborate  it.  The  remaining 
records  were  then  combined  and  Idealized  as  follows. 

Data  later  than  about  three  microseconds  after  shock  arrival  were 
discarded  because  of  the  cross  talk  from  downstream  gages  discussed  In 
the  previous  subsection.  (The  profiles  were  extrapolated  slightly 
beyond  the  first  occurrence  of  the  perturbations  if  we  judged  the  data 
warranted  It.)  Redundant  profiles  were  then  averaged  and  smoothed  by 
hand.  Finally,  the  four  resulting  histories  were  Idealized  In  two 
ways.  In  the  first,  the  dip  occurring  after  the  first  rapid  spike  was 
Included;  In  the  second.  It  was  omitted.  (The  dip  was  omitted  by 
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extrapolating  the  late  time,  slow  decay  portion  of  the  history  backward 
In  time  until  the  falling  edge  of  the  spike  was  intersected.)  Although 
the  dip  is  quite  apparent  In  the  original  oscillographs  (Figure  8)  and 
the  digitized  histories  (Figures  9  and  10),  we  chose  to  omit  it  in  one 
idealization  because  we  were  not  sure  we  would  be  able  to  incorporate  it 
in  a  reactive  hydrodynamic  model  of  the  detonation  process  in  NA/NP. 

The  resulting  profiles  for  both  idealizations  are  shown  in  Figure  13. 

We  were  now  in  a  position  to  characterize  the  reactive  flow 

observed  in  NA/NP.  Both  idealizations  included  a  spike  propagating  at 

the  steady  velocity  of  6.48  km/s  followed  by  a  slightly  dispersive 

rarefaction,  in  conformity  with  the  qualitative  behavior  expected  based 

on  the  ZND  detonation  model.  The  detonation  velocity  of  6.48  km/s  is  in 

excellent  agreement  with  the  value  of  6.50  km/s  inferred  from  a  plot  of 
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detonation  velocity  versus  acid  concentration  presented  by  Mallory. 

The  precise  starting  point  for  the  nonsteady  flow,  that  is,  the  CJ 
point,  cannot  be  resolved.  Therefore,  for  the  idealization  including 
the  dip,  we  arbitrarily  select  it  to  be  the  maximum  particle  velocity 
reached  after  the  dip,  and  for  the  idealization  omitting  the  dip,  we 
choose  it  to  be  the  apparent  discontinuity  in  decay  rate  after  the 
spike.  A  surprising  result  was  that  NA/NP  has  a  significant  reaction 
time:  0.2  ps  for  idealization  omitting  the  dip  and  0.55  ns  with  the  dip 
included  and  interpreted  as  being  in  the  reaction  zone. 

The  two  idealized  data  sets  shown  in  Figure  13  were  Lagrange 
analyzed  using  the  SRI  GUINSY  code  discussed  in  Section  2.2  The  results 
are  shown  in  the  p-v  plane  in  Figure  14.  This  figure  represents  the 
loci  of  p-v  paths  followed  by  four  different  particles  having  the 
Lagrange  coordinates  of  the  last  four  gage  planes  in  the  second  NA/NP 
experiment.  The  CJ  point  indicated  is  an  average  of  the  values  on  the 
middle  two  gage  planes  for  the  two  idealizations. 

Figure  14  represents  the  results  of  the  present  Lagrange  gage  and 
analysis  characterization  of  the  steady  detonation  process  in  NA/NP.  It 
shows  the  Rayleigh  line  followed  in  the  shock  compression  to  the  maximum 
pressure  resolved  by  the  gages  (less  than  the  ZND  spike  pressure  because 
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Figure  14.  Lagrange  analysis  calculations  of  pressure-specific 
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of  the  response  time  of  our  EPV  gages),  the  reactive  release  states  also 
lying  on  the  Rayleigh  line,  the  CJ  state,  and  the  initial  adiabatic 
release  states.  The  CJ  state  variables  are  D  =  6.48  km/s,  p  =  13.3  GPa, 

V  =  0.5986  m^/Mg,  V  =  0.7494,  u  =  1.64  km/s,  and  the  adiabatic  exponent 
(logarithmic  derivative  of  isentrope)  y  =  2.99. 

3.4  TIGER  CALCULATIONS 

TIGER  calculations  using  the  BKW  and  the  JCZ3  EOSs  were  performed 
to  extend  the  NA/NP  CJ  adiabatic  release  path  to  large  volumes  by  model¬ 
ing  the  release  process  as  being  isentropic.  First,  we  used  both  EOSs 
to  predict  the  CJ  state  as  well  as  the  CT  release  isentrope.  The 
results  in  the  p-v  plane  are  shown  in  Figure  15,  on  a  log-log  scale, 
along  with  the  Lagrange  gage  and  analysis  CJ  state  and  a  release  path 
discussed  later  in  this  section. 

The  CJ  state  variables  from  the  two  calculations  are  listed 
below.  From  the  BKW  calculation,  D  =  7.14  km/s,  p  =  16.40  GPa, 

V  =  0.5937  m^/Mg,  u  =  1.83  km/s,  and  y  =  2.90.  From  the  JCZ3  calcu¬ 
lations,  D  =  6.12  km/s,  p  =  11.30  GPa,  v  =  0.6060  m^/Mg,  u  =  1.48  km/s, 
and  Y  “  3.14. 

The  CJ  state  calculated  by  JCZ3  more  closely  matches  the  Lagrange 
gage  and  analysis  results.  A  more  significant  difference,  however,  for 
our  purposes,  is  that  the  JCZ3  calculation  extends  to  large  volume 
expansions  so  we  can  apply  the  JWL  fitting  procedures  described  in 
Section  2,  whereas  the  BKW  calculation  does  not.  The  difference  is 
believed  to  result  from  the  construction  of  TIGER  and  its  interaction 
with  this  particular  detonation  problem  and  not  to  reflect  any  diffi¬ 
culty  with  the  BKW  EOS.  Because  of  the  limited  range  of  the  present  BKW 
calculations,  we  elected  to  use  JCZ3,  in  the  manner  described  below,  to 
match  the  Lagrange  gage  and  analysis  results  and  extend  the  CJ  release 
adiabat  to  large  volumes. 

We  matched  the  experimental  results  simply  by  initiating  the 
TIGER/JCZ3  Isentrope  calculations  at  a  p-v  point  on  the  experimental 
release  path.  We  used  both  the  CJ  point  and  the  final  point  on  the  CJ 
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release  adlabat.  In  the  regions  of  overlap  there  was  no  detectable 
difference  between  either  isentrope  or  the  experimental  release  path. 
Both  isentropes  are  shown  in  Figure  15  and  are  indistinguishable.  For 
computational  convenience,  we  selected  the  Isentrope  centered  on  the 
experimental  CJ  state  to  generate  initial  JWL  coefficients  for  NA/NP. 

3.5  CALCULATION  OF  JWL  COEFFICIENTS 

Initial  JWL  coefficients  for  NA/NP  were  calculated  by  applying  the 
method  described  in  Section  2  to  the  basis  NA/NP  CJ  release  adlabat 
taken  to  be  the  JCZ3  CJ  Isentrope  through  the  experimental  CJ  state 
(Figure  15).  The  value  was  30.  The  resulting  JWL  coefficients  are 
listed  in  Table  1. 


Table  1.  NA/NP  JWL  coefficients. 


Initial 

State 

CJ  State 

Equatlon-of-State  Paraaeters 

Po 

P 

D 

Eo 

r 

A 

B 

C 

Ri 

R2 

(0 

(Mg/a^) 

(Mbar) 

(cm/iis) 

(Mbar) 

(Mbar) 

(Mbar) 

(Mbar) 

1.252 

0.133 

0.648 

0.0697 

2.991 

7.554 

0.1962 

0.012717 

6.366 

1.810 

0.285 

The  JWL  and  basis  CJ  release  paths  in  the  p-V  plane  are  shown  in 
Figure  16.  Figure  17  shows  the  adiabatic  expansion  coefficients,  y»  the 
logarithmic  derivatives  of  the  release  paths  in  Figure  16.  The  cor¬ 
relation  between  the  Isentropes  in  Figure  16  is  controlled  by  the  fit¬ 
ting  process;  the  behavior  in  the  y-V  plane.  Figure  17,  is  not  con¬ 
strained  . 

The  NA/NP  JWL  coefficients  in  Table  1  are  the  major  result  of  the 
research  described  in  this  section  and  can  be  used  in  DNA  hydrocode 
calculations  to  predict  the  performance  of  NA/NP  in  NWE  simulators. 
Additional  discussion  is  given  in  the  next  subsection,  and  recommenda¬ 
tions  regarding  the  use  of  NA/NP  and  further  research  are  given  in 
Section  5. 
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Figure  16.  JWL  CJ  release  isentrope  for  NA/NP. 


Table  2.  Comparison  of  NA/NP  and  NM  JWL  coefficients. 

MA/NP 


Initial 

State 

CJ  State 

Equatlon-of-State  Parameters 

PQ 

(Mg/m^) 

P 

(Mbar) 

0 

(cm/^s) 

Eo 

(Mbar) 

■ 

A 

(Mbar) 

B 

(Mbar) 

C 

(Mbar) 

h 

R2 

(a) 

1.252 

0.133 

0.6A8 

0.0697 

2.991 

7.554 

0.1962 

0.012717 

6.366 

1.810 

0.285 

1.128 

0.125 

0.628 

0.0510 

2.538 

2.092 

0.05689 

0.00770 

4.40 

1.20 

0.30 

The  JWL  coefficients  generated  here  for  NA/NP  differ  from  typical 
cylinder  test  values  in  that  the  nonlinear  coefficients,  and  R2»  are 
larger.  We  believe  this  is  because  of  the  larger  expansion  range  over 
which  the  Isentrope  is  constrained  here  and  possibly  because  of  the 
I  difference  in  the  energy  constraints  between  the  two  methods.  The 

present  coefficients  are  within  the  range  of  results  from  cylinder  tests 
and  are  believed  to  be  appropriate  for  a  range  of  DNA  applications  as 
Intended.  For  specific  applications,  specialized  coefficients  could  be 
generated  from  dedicated  experiments  simulating  the  application,  if 
necessary. 

I 

'  The  Lagrange  gage  and  analysis  results  indicated  a  surprisingly 

long  reaction  zone,  between  0.2  and  0.6  ps,  apparently  with  resolvable 
structure,  for  NA/NP  undergoing  steady  planar  detonation.  This  may 
I  result  from  the  effects  of  the  water  in  the  reactants,  which  may  Impede 

the  kinetics  of  the  chemical  reaction  and  will  also  undergo  a  phase 

I 

j  change  and  affect  the  temperature  distribution  in  the  products. 

The  SDT  experiment  was  ambiguous.  The  results  were  not  repro- 

[  ducible  in  that  gage  records  1  through  5  differ  from  gage  records  6 

through  10.  Also,  over  the  range  spanned  by  the  gages  (50  to  100  mm 

I  from  the  initiator  system)  the  input  shock  had  neither  built  to  steady 

J  detonation  nor  decayed.  Some  of  the  particle  velocity  histories  showed 

• 

«  a  late  time  pulse  building  with  depth  in  the  target,  behind  a  narrow 

»  - - - 

From  Reference  12. 

! 
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relatively  constant  Initial  spike.  This  structure  has  a  qualitative 
resemblance  to  the  profiles  observed  in  the  steady  detonation  experi¬ 
ment.  The  results  indicate  that  the  shock  initiation  of  detonation  in 
NA/NP  may  be  unusual,  and  in  some  cases  sluggish,  and  should  be  Investl 
gated  further. 

We  believe  that  NA/NP  is  a  good  liquid  HE  offering  an  attractive 
low  cost  alternative  to  NM.  Neither  the  Baume  42  weak  nitric  acid  nor 
the  nitropropane  presented  undue  handling  problems. 


SECTION  4 


NPN  CHARACTERIZATION 


Initial  JViL  equatlon-of-state  parameters  were  developed  for  the 
detonation  products  of  NPN,  a  blasting  agent  somewhat  like  ANFO.  The 
specific  composition  and  Initial  density  of  the  NPN  selected  for  charac¬ 
terization  were  chosen  to  be  similar  to  those  for  the  NPN  Investigated 
previously  for  DNA  in  an  alrblast  performance  test  by  Franzen  and 
Wisotski. 

The  tasks  performed  included  selection  and  fabrication  of  a  parti¬ 
cular  NPN  study  material;  Lagrange  gage  experiments  and  analysis  to 
determine  shock  compression,  CJ,  and  Initial  CJ  adiabatic  release 
states;  TIGER  calculations  to  interpret  the  experiments  and  extend  the 
experimental  CJ  release  adlabat  to  large  volume  expansions;  and  deter¬ 
mination  of  JVL  coefficients.  The  results  of  these  tasks  are  described 
In  the  following  subsections. 

4.1  DESCRIPTION  OF  MATERIAL 
Background 

NPN  (nltropropane  nitrate)  is  composed  of  the  solid  oxidizer 

ammonium  nitrate  (AN)  and  the  liquid  fuels  nltropropane  (NP)  and  methyl 

alcohol  (MA).  The  material  was  developed  initially  by  John  R.  Post  as 

an  alternative  blasting  agent  to  ANFO  with  improved  performance  for 

2 1 

mining  applications.  Basic  NPN  physical  characteristics  and  initial 

performance  data,  as  well  as  a  fairly  standardized  baseline  composition, 

22 

were  described  by  Post  In  1977  and  are  summarized  below. 

The  baseline  NPN  composition  given  by  Post  was  87  wt%  AN  and  13  wt% 
fuel  (NP  and  MA  In  equal  volumes).  This  composition  was  reported  to  be 
oxygen-balanced  and  could  be  loaded  to  densities  ranging  from  about  1  to 
1.4  Mg/ra  depending  on  the  degree  of  compaction  induced  by  the  loading 
procedures  and  by  its  own  weight.  Approximate  detonation  velocities 


observed  by  Post  were  between  4.5  and  6  km/s  depending  on  Initial  den¬ 
sity  and  confinement. 

We  are  aware  of  only  two  experimental  investigations  of  NPN  other 
than  Post's  work  cited  above.  In  1977  Mallory  measured  the  detonation 
velocity  in  thick-wall  1-lnch-ID  steel  pipes  for  the  baseline  composi¬ 
tion  at  an  initial  density  of  1.298  Mg/m^.^^  He  observed  a  detonation 
velocity  of  5.206  ±  0.028  km/s.  In  1983  Franzen  and  Wlsotskl^*^  per¬ 
formed  a  100-pound  quasl-spherlcal  airblast  test  for  an  NPN  composition 
slightly  modified  from  the  baseline  material  (see  next  subsection)  at  an 
initial  density  of  1.05  Mg/m^.  They  observed  propagation  velocities 
between  2.8  and  4.1  km/s  in  the  detonation  buildup  regime. 

In  addition  to  the  specific  investigations  cited  above,  two  recent 
DNA  reports  review  useful  NPN  information.  The  first  is  an  informal 

O  / 

report  prepared  by  Franzen  summarizing  available  documents  on  NPN.^ 

The  other  is  a  technical  report  by  McMullan^^  containing  an  evaluation 
of  the  properties  of  a  number  of  explosives,  including  those  of  NPN 
summarized  by  Franzen^^,  for  large-scale  ONA  air  blast  tests. 

Present  Study  Material 

The  composition  of  the  present  study  material,  like  the  material 
studied  by  Franzen  and  Wlsotski,  is  84.8  wt%  AN,  0.2  wt%  raethocel  (a 
stabilizer  to  prevent  migration  or  puddling  of  the  liquid  fuels, and 
15.0  wt%  fuel.  However,  the  present  study  material  differs  slightly 
from  that  investigated  by  Franzen  and  Wlsotski  in  that  the  relative 
proportions  of  the  fuel  components,  MA  and  NP,  are  equal  by  weight,  not 
volume.  Both  compositions,  which  can  be  considered  equivalent  for  most 
purposes,  have  a  slightly  higher  fuel  concentration  than  the  oxygen- 
balanced  baseline  material.  Post^^  recommended  this  overfueling  for 
both  this  study  and  the  previous  study  by  Franzen  and  Wlsotski  to 
Increase  the  dampness-induced  cohesiveness  and  therefore  the  achievable 
initial  density  in  the  NPN.  In  both  studies,  the  fragility  of  the 
embedded  sensors  precluded  the  use  of  heavy  tamping  to  achieve  the  high 
NPN  loading  density  characteristic  of  field  applications. 
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The  NPN  tested  in  this  work  was  constructed  from  the  following 
materials.  The  AN  was  crushed,  porous,  uncoated,  fertilizer-grade 
prills.  Crushing  the  rather  uniform  spherical  AN  prills  increases  both 
packing  densities  and  intimacy  of  mixing  of  the  AN  and  the  fuels. 

Because  degree  of  crushing  and  prill  porosity  are  AN  parameters  that  may 
affect  performance  and/or  sensitivity  of  NPN,  we  purchased  the  crushed 
AN  (premixed  with  methocel)  directly  from  Post,  who  also  supplied  these 
materials  for  the  work  by  Franzen  and  Wlsotskl.  The  approximate  screen 
analysis  provided  by  Post  for  our  crushed  AN,  and  described  as  represen¬ 
tative  of  that  obtained  for  AN  crushed  in  the  field  with  commercial 
hammer  mills,  is  tabulated  below.  ® 


Mesh  wt% 


+20 

30.5 

-20 

+60 

40.5 

-60 

+100 

11.8 

•100 

17.2 

The  NP  was  furnished  by  International  Minerals  and  Chemical  Corporation, 
IMC  Chemical  Group,  and  is  a  mixture  of  the  two  stable  Isomers.  The  MA 
was  standard  reagent  grade  material. 

Procedure  for  Producing  Uniform  Density  Specimens 

Before  the  Lagrange  particle  velocity  gage  experiments  were  per¬ 
formed,  a  significant  effort  was  made  to  identify  a  method  of  loading 
the  NPN  specimen  material  into  the  Lagrange  gage  target  chamber  so  that: 

(1)  The  NPN,  which  is  highly  compactible,  would  have  a 
constant  density  throughout  the  chamber,  particularly 
in  the  vicinity  of  the  gages. 

(2)  The  density  would  be  close  to  1.05  Mg/m^. 

(3)  The  fragile  particle  velocity  gages  would  not  be 
distorted  and  would  be  in  intimate  contact  with  the 
NPN. 


To  this  end,  a  "sampling  densitometer”  was  constructed  in  the  shape 
of  a  cylindrical  chamber  with  an  inside  diameter  of  0.124  m  and  a  height 


of  approximately  0.30  m.  The  densitometer  Incorporated  thin  steel 
sheets  positioned  at  four  heights  in  the  cylinder  that  could  be  pulled 
across  guillotine-style  to  physically  separate  the  NPN  into  a  set  of 
shorter  cylindrical  chambers  with  minimal  sample  disturbance.  These 
sections  could  be  weighed  separately  after  various  trial  loading  proce¬ 
dures  to  obtain  the  average  NPN  densities  at  heights  within  the  densito¬ 
meter  corresponding  to  the  heights  of  the  gage  elements  in  the  Lagrange 
gage  test  chamber. 

Figure  19  shows  the  densitometer  results  for  the  five  filling 
methods  methods  that  were  evaluated: 

(1)  Strong  vertical  shaking  of  the  densitometer  during 
and  after  loading  the  NPN  through  a  large  funnel. 

(2)  Pneumatically  vibrating  the  chamber  during  and  after 
loading  the  NPN  through  a  large  funnel. 

(3)  Drop  loading  the  NPN  through  a  vibrating  mesh  screen 
located  about  one  meter  above  the  chamber. 

(4)  Same  as  method  (3),  followed  by  vertically  compres¬ 
sing  the  sample  approximately  2.5  cm. 

(5)  Same  as  method  (3),  followed  by  deadweight  compres¬ 
sing  the  sample  to  a  predetermined  density,  while 
pneumatically  vibrating  the  chamber. 

As  seen  from  Figure  19,  the  first  two  methods  resulted  In  bulk 
densities  averaged  over  the  entire  densitometer  that  were  reasonably 
close  to  1.05  Mg/m^,  but  with  large  density  gradients  from  top  to  bot¬ 
tom.  The  last  three  methods  resulted  In  lower  average  densities,  pre¬ 
sumably  because  the  screen  mesh  reduced  the  Initial  compacting  ten¬ 
dencies  of  the  material.  Method  (5),  In  which  the  deadweight  compres¬ 
sion  during  vibration  reduced  the  pressure  gradient  in  the  sample  while 
allowing  the  material  to  flow,  resulted  In  a  very  low  density  gradient, 
at  an  average  density  of  1.01  Mg/m^.  Attempts  to  Increase  the  bulk 
density  while  retaining  a  low  density  gradient  were  unsuccessful. 

Method  (5)  was  chosen  as  the  basis  for  filling  the  chambers  for  the 
particle  velocity  gage  tests. 


HEIGHT  (cm) 


O  Method  1 
•  Method  2 
A  Method  3 


▲  Method  4 
■  Method  5 


0.94  0.98  1.02  1.06  1.10 
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Region  of  Particle 
Velocity  Gage 
Elements  in  the 
Lagrangian  Gage 
Test  Chamber 
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Figure  19.  Density  of  NPN  at  various  heights  in  the  densitometer  chamber 
for  five  filling  procedures. 
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4.2  LAGRANGE  GAGE  EXPERIMENTS 


Lagrange  Gage  Target  Chambers  and  Filling  Procedures 

The  Lagrange  gage  target  chambers.  Figure  20,  were  PVC  cylinders 
approximately  0.30  m  OD,  0.28  m  ID,  and  0.15  m  high.  The  chambers  were 
sealed  at  the  bottom  with  a  PMMA  base  plate  containing  a  standard  par¬ 
ticle  velocity  gage  block  of  the  type  described  in  Section  2  (see 
Figures  3  and  4)  and  used  in  the  NA/NP  experiments  described  in  Section 
3.  An  auxiliary  0.105-m-dlameter  filling  hole  with  a  cap  (for  use  in 
the  filling  operation  described  below)  was  cut  in  the  cylindrical  wall 
of  each  PVC  chamber  in  line  with  the  axis  of  the  gage  block. 

The  gage  blocks,  each  containing  the  usual  five  pairs  of  0.15-mm- 
thlck  aluminum  foil  particle  velocity  gages,  were  mounted  in  the  base 
plates  so  that  the  active  elements  of  the  gages  were  parallel  to  the  top 
surface  of  the  chambers  to  within  ±0.025  mm.  The  active  elements  of  the 
tallest  gage  pair  (that  is,  the  first  pair  to  enter  the  flow)  were 
located  82.5  mm  below  the  upper  NPN  surface,  and  the  active  elements  of 
the  shortest  gage  pair  (the  last  pair  to  enter  the  flow)  were  50.8  mm 
below  that.  (The  difference  in  heights  between  successive  gage  pairs  is 
12.7  mm  as  discussed  in  Section  2.)  T’  .•  inner  surface  of  the  base 
plates  and  the  tops  of  the  gage  blocks  were  coated  with  paraffin  for 
protection  from  the  solvent  action  of  the  NP.  Each  target  chamber  was 
loaded  immediately  before  the  experiment  and  covered  with  a  1.7-imn-thlck 
PMMA  plate  to  prevent  evaporative  loss  of  the  fuel  and  to  separate  the 
driver  charge  from  the  NPN  solvents. 

The  target  chamber  filling  procedure  was  based  on  the  principle  of 
method  (5)  developed  with  the  sampling  densitometer,  modified  to 
minimize  disturbance  to  the  particle  velocity  gages.  The  preweighed 
components  were  mixed  to  produce  a  slightly  damp  compressible  powder  as 
before.  Approximately  90%  of  the  material  required  to  achieve  the 
desired  density  was  drop  loaded  into  the  pneumatically  vibrated  test 
chamber  from  the  top  through  a  mesh  screen.  The  top  of  the  chamber  was 
then  sealed  with  the  PMMA  cover  plate  and  a  temporary  PMMA  backing 
plate.  The  chamber  was  rotated  90°  until  the  long  dimension  of  the  gage 
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Figure  20.  NPN  Lagrange  gage  target  chamber  (PVC) 

(Top  and  bottom  are  sealed  with  PMMA  plates.  The  four 
blocks  beneath  chamber  are  temporary  supports  for  this 
photograph. 
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block  was  vertical,  and  the  remaining  material  was  then  deadweight  com¬ 
pressed  through  the  auxiliary  fill  hole  as  in  method  (5).  The  motion  of 
the  dense  NPN  in  the  final  stages  of  filling  was  thus  In  the  direction 
of  maximum  gage  element  rigidity  (parallel  to  the  plane  of  the  gage 
foil). 

For  the  first  Lagrange  gage  experiment,  the  chamber  was  subjected 
to  planetary  pneumatic  vibration  while  the  main  portion  of  the  NPN  was 
added  and  before  the  chamber  was  rotated  90°.  After  the  experiment  we 
discovered  that  the  planetary  vibration  caused  bulk  rotation  of  the  NPN 
within  the  circular  cross-section  chamber  and  severely  distorted  the 
gages.  For  the  second  experiment,  linear  vibration  was  used  before  the 
chamber  was  tilted,  and  no  gage  distortion  was  detected  or  is  believed 
to  have  occurred. 

The  NPN  specimen  uniformity  In  the  Lagrange  gage  experiments  Is 
estimated  from  the  results  of  the  sampling  densitometer  study,  from 
visual  Inspection  of  the  specimen  material  through  the  transparent 
target  chamber  cover  during  and  after  loading,  from  experience  gained  in 
trial  loadings  of  the  target  chambers,  and  from  the  results  of  the 
Lagrange  gage  experiments  themselves.  The  final  average  (bulk)  density 
was  1.01  Mg/ra^,  and  we  estimate  that  the  density  variation  throughout 
the  majority  of  the  chamber  was  less  than  ±2%.  The  NPN  appeared  to  be 
In  Intimate  contact  with  the  gages. 

During  the  final  step  of  the  chamber  loading  operation  for  the 
second  experiment,  we  noted  that  small  gaps  developed  along  the  upper 
portion  of  the  cylindrical  Interface  between  the  NPN  and  the  target 
chamber  wall.  Immediately  adjacent  to  the  filling  tube.  These  resulted 
from  the  combined  effects  of  settling  and  nonhydrodynamlc  rheology  in 
the  NPN.  We  believe  these  gaps  had  a  very  small  effect  on  the  bulk 
properties  of  the  NPN;  however,  they  did  asyraetrlcally  alter  the  target 
width  such  that  the  effective  target  radius  was  smaller  on  the  side  near 
the  filling  hole  than  on  the  other  side.  In  spite  of  this  difficulty, 
we  believe  that  the  specimen  density  in  the  Lagrange  gage  experiments  is 
probably  more  uniform  than  could  be  achieved  in  the  field. 


Initiation  System 

We  selected  Baratol  as  an  appropriate  driver  for  rapidly  Initiating 
steady  detonation  In  NPN.  Slight  underdrlvlng  was  assumed  to  be  the 
optimal  condition  for  Initiating  steady  detonation  In  the  present 
Lagrange  gage  experiments.  We  assumed  that  the  reaction  time  In  any 
candidate  military  HE  driver  was  negligible  compared  with  that  In  the 
NPN.  We  then  sought  a  driver  for  which  the  shock  Impedance  mismatch 
between  the  driver  detonation  products  and  the  unreacted  NPN  produced  a 
pressure  In  the  unreacted  NPN  about  equal  to  the  expected  NPN  CJ  pres¬ 
sure.  Because  this  pressure  on  the  unreacted  NPN  Hugonlot  Is  well  below 
the  expected  spike  pressure  on  the  Hugonlot  during  steady  detonation  and 
because  of  the  long  postulated  NPN  steady  reaction  time,  this  Initiation 
pulse  Is  taken  to  underdrive  NPN  detonation.  However,  because  of  the 
long  duration  of  the  driver  Induced  pulse  and  its  significant  percentage 
of  the  expected  NPN  steady  detonation  spike  amplitude,  this  Initiation 
pulse  was  estimated  to  satisfy  any  threshhold  criteria  for  shock  Initia¬ 
tion  of  detonation  and  thus  to  provide  the  slight  underdrlvlng  but 
reliable  Initiation  we  sought. 

We  estimated  by  an  Impedance  mismatch  calculation  that  Baratol 
would  generate  an  Input  shock  pressure  In  unreacted  NPN  of  about  7.0*1 
GPa.  For  this  estimate  we  very  roughly  approximated  the  Hugonlot  for 
unreacted  NPN  from  unreacted  AN  data.^^  Since  7  GPa  was  the  CJ  pressure 
we  estimated  for  NPN  at  a  density  of  about  1  Mg/m^  from  data  by  Post,  we 
selected  Baratol  for  the  HE  driver. 


Experimental  Results 

Two  NPN  Lagrange  gage  detonation  experiments  were  performed.  The 
purpose  of  the  first  was  to  evaluate  gage/specimen  compatibility  and  to 
determine  optimal  Instrumentation  settings.  In  this  experiment  a  20-cm- 
dlameter  P-80  plane  wave  lens  was  used  to  Initiate  detonation.  Since 
the  base  charge  In  our  P-80  lenses  Is  Baratol,  the  amplitude  of  the 
shock  Initiation  pulse  driven  Into  the  NPN  target  Is  that  resulting  from 
the  Impedance  mismatch  between  Baratol  detonation  products  (considering 
the  reaction  zone  In  Baratol  to  be  negligible)  and  unreacted  NPN. 


Data  were  obtained  from  all  gages  In  the  first  experiment,  result¬ 
ing  In  good  time-of-arrlval  Information  for  optimizing  oscilloscope 
sweep  speeds  and  delays  in  the  second  experiment.  However,  the  shapes 
of  the  recorded  profiles  were  Irregular  and  Inconsistent,  preventing 
quantitative  analysis  of  the  detonation  process  in  NPN.  As  discussed 
previously,  these  anomalous  profiles  are  attributed  to  the  gage  distor¬ 
tions  resulting  from  the  planetary  vibration  used  in  filling  the  target 
chamber  for  this  experiment.  Modified  filling  procedures  were  used  in 
the  second  experiment. 

The  purpose  of  the  second  Lagrange  gage  experiment  was  to  provide 
data  on  the  detonation  process  in  NPN.  In  this  experiment,  a  30-cm- 
diameter  plane  wave  lens  (P-120)  and  a  10-cm-thlck  by  30-cm-d lame ter 
Baratol  pad  were  used  to  initiate  the  NPN  to  increase  the  uniaxial 
strain  recording  time  at  the  gage  locations. 

Records  were  obtained  from  all  particle  velocity  gages.  The  oscil¬ 
lographs  are  shown  in  Figure  21.  For  these  oscillographs,  vertical 
sensitivities  are  500  mV/divlslon  and  horizontal  sensitivities  are 
1  ps/dlvlslon.  Short  blanks  on  all  the  scope  records  are  simultaneous 
fiducial  marks  generated  by  a  repetitive  "z-axls  blanking"  beam  inten¬ 
sity  modulation  system  used  to  time-correlate  the  oscillographs. 

The  NPN  gage  records  in  Figure  21  are  judged  to  be  better  than 
average  detonation  records  in  two  ways  but  worse  in  a  third.  We  have 
previously  used  this  gage  design  to  observe  detonation  in  higher  perfor- 
mance  explosives  Including  Composition  B,  Amatex  20,*°  and  the 
eutectic  explosive  EAK,^^  as  well  as  in  the  NA/NP  experiments  described 
in  Section  3.  Typical  results.  Composition  B  records  from  reference  27, 
are  reproduced  in  Figure  22.  The  records  in  Figure  21  show  considerably 
less  high  frequency  noise  on  the  first  gage  plane  than  the  records  from 
the  other  materials  such  as  the  Composition  B  data  in  Figure  22  (the 
NA/NP  results.  Figure  8,  are  also  free  of  this  noise).  In  addition,  the 
records  in  Figure  21  contain  considerably  smaller  perturbations  on 
upstream  records  when  downstream  gages  enter  the  shock  front  than  we  see 
in  the  other  materials.  These  two  Improvements  may  result  from  the 
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lower  temperatures,  pressures,  and  conductivities  expected  in  the  NPN 
reaction  zone  than  in  the  reaction  zones  of  the  higher  performance 
materials. 

However,  the  NPN  records  in  Figure  21  contain  more  pronounced 
Intermediate  frequency,  low  amplitude  oscillations  during  stress  release 
than  the  records  from  the  other  materials  in  this  region.  The  origin  of 
this  structure  is  unknown,  but  it  may  be  related  either  to  the  particu¬ 
late  nature  of  the  NPN  or  to  the  relatively  large  plastic  strains 
experienced  by  the  leads  in  the  present  experiments. 

The  ten  particle  velocity  gage  oscilloscope  records  in  Figure  21 
were  digitized  and  reduced  to  the  particle  velocity  histories  shown 
together  on  a  common  time  scale  in  Figure  23.  This  figure  shows  the 
overall  nature  of  the  flow  in  a  single  composite  plot.  In  the  following 
section  these  histories  evaluated  in  various  subgroups  in  which  their 
individual  features  are  more  readily  viewed. 

Recall  that  the  gages  are  numbered  sequentially  from  one  side  of 
the  gage  block  to  the  other  (see  Figure  4).  Thus,  each  set  of  histories 
(1,10),  (2,9),  (3,8),  (4,7),  and  (5,6)  results  from  a  pair  of  redundant 
particle  velocity  measurements  at  one  Lagrange  coordinate,  and  the  five 
Lagrange  coordinates  at  which  the  measurements  are  made  are  initially 
12.7  mm  apart  in  the  wave  propagation  direction.  The  gage  pair  (1,10) 
is  the  closest  to  the  driver  and  thus  the  first  to  enter  the  flow;  the 
pair  (5,6)  is  the  closest  to  the  rear  boundary  of  the  target  and  thus 
the  last  to  enter  the  flow. 

4.3  DATA  ANALYSIS 

The  ten  particle  velocity  records  in  Figure  23  were  first  compared 
to  (1)  obtain  a  best  nonredundant  set  of  histories  for  Lagrange  analysis 
and  (2)  characterize  the  reactive  flow  in  NPN.  We  performed  three  types 
of  comparisons:  Intra-slde,  intra-pair,  and  inter-side.  In  the  intra- 
side  comparisons,  we  considered  the  five  profiles  from  each  side  of  the 
target  as  independent  groups  to  evaluate  the  data  and  the  flow  in  each 
half  of  the  target  as  separate  experiments.  In  the  Intra-palr 
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comparisons,  we  compared  the  two  redundant  measurements  at  each  of  the 
five  Lagrange  coordinates  to  determine  measurement  reproducibility  at 
these  locations*  In  the  Inter-side  comparison,  we  compared  the  five 
records  from  one  side  of  the  target  with  those  from  the  other  to 
evaluate  overall  target  response.  Below  we  summarize  the  comparisons, 
select  a  data  set  for  analysis,  and  give  an  Initial  characterization  of 
the  flow.  Additional  interpretation  is  given  after  the  thermodynamic 
TIGER  calculations  are  presented. 

Intra-side  Comparisons 

Histories  1-5,  recorded  by  the  gages  that  were  on  the  top  side  of 
the  target  chamber  during  the  final  filling  operation,  are  shown  in 
Figure  24.  The  peak  particle  velocity  amplitudes  increase  monotonically 
with  propagation  distance  in  the  region  spanned  by  the  first  three  gages 
(1-3)  and  remain  constant,  within  experimental  error,  in  the  region 
spanned  by  the  last  three  (3-5).  This  is  consistent  with  steady  flow 
beginning  near  gage  3. 

The  shapes  of  histories  3-5  are  essentially  identical  within 
experimental  error,  again  consistent  with  steady  flow  in  the  region 
spanned  by  these  gages.  (Experimental  accuracy  of  the  particle  velocity 
measurements  is  estimated  to  be  i5%.)  History  1  has  an  anomalous  shape 
compared  with  histories  3-5  and,  for  the  first  microsecond  after  enter¬ 
ing  the  flow,  is  almost  50%  lower  in  amplitude  than  these  histories. 
History  2  is  for  the  most  part  similar  in  shape  to  3-5,  but  is  as  much 
as  20%  lower  in  amplitude  during  the  first  microsecond. 

The  histories  from  the  bottom  half  of  the  target,  6-10,  are  shown 
in  Figure  25.  The  peak  amplitudes  and  shapes  of  the  histories  from  the 
last  four  gages  (9-6)  are  consistent  with  steady  flow  in  the  region 
spanned  by  these  gages.  No  growth  of  peak  amplitude  is  observed  on  the 
first  gages;  in  fact,  the  first  profile,  10,  is  anomalously  large  at 
early  times. 


The  shock  transit  velocities  between  adjacent  gages  (not  redundant 
gage  pairs),  for  each  side  of  the  target,  and  the  mean  transit 


velocities  and  standard  deviations  for  each  side,  are  listed  in  Table  3 
along  with  the  shock  arrival  times  from  which  they  were  calculated. 

Here  arrival  times  are  for  the  midpoint  of  the  leading  edge  of  the 
compressive  wave,  following  recommended  shock  wave  practice.  In  the 
side  of  the  target  away  from  the  auxiliary  filler,  the  transit  veloc¬ 
ities  agree  within  measurement  accuracy  (estimated  to  be  2%-3%).  This 
is  consistent  with  steady  flow  propagating  at  5.07  km/s.  (We  note  here, 
however,  that  for  reasons  presented  in  the  following  discussion,  the 
steady  velocity  attributed  to  the  detonation  wave  in  NPN  is  5.08 
km/s.)  In  the  side  of  the  target  nearer  the  auxiliary  filling  hole,  the 
shock  propagation  velocity  differences  exceed  experimental  error, 
indicating  less  uniform  flow. 

Intra-pair  Comparisons 

Figure  26  shows  the  redundant  pairs  of  histories  measured  at  each 
of  the  five  Lagrange  gage  coordinates.  The  two  records  from  the  first 
gage  pair  to  enter  the  flow,  1  and  10,  are  qualitatively  different. 

Note  that  history  1  is  both  lower  in  initial  particle  velocity  amplitude 
and  later  in  time  of  first  motion  than  history  10. 


The  records  in  the  remaining  pairs  are  qualitatively  similar  to 
each  other.  However,  in  the  second  pair,  2  and  9,  2  is  initially  15%- 
20%  lower  than  9  and  later  is  higher  by  a  similar  percentage.  These 
differences  exceed  the  estimated  experimental  errors.  In  the  third 
pair,  3  and  8,  the  differences  during  the  shock  compression  also  exceed 
the  estimated  experimental  errors,  but  are  the  result  of  a  time  correla¬ 
tion  difference  rather  than  a  profile  shape  difference.  If  3  is  shifted 
back  approximately  0.1  microsecond,  reducing  the  scatter  in  shock  propa¬ 
gation  velocity,  these  records  agree  to  within  experimental  error  every¬ 
where.  The  records  in  each  of  the  last  two  pairs  to  enter  the  flow 
generally  agree  to  much  better  than  ±5%,  indicating  that  both  the  flow 
and  the  measurements  are  reproducible  at  these  locations. 


Table  3. 


Shock  arrival  times,  transit  velocities,  and  detonation 
velocity  for  second  NPN  experiment® 


Shock  Arrival 

Times  at  Gages 

Gage  Pair  Transit  Velocities 

Gage 

Shock 

Arrival 

Time'’ 

(l^s) 

Side  Near  Filler  | 

Side  Away  From  Filler  | 

Gage 

Pair 

Transit 

Velocity 

(km/s) 

Gage 

Pair 

Transit 

Velocity 

(km/s) 

1 

(1.23) 

10 

(0.90) 

1-2 

(5.75) 

10-9 

(5.04) 

2 

(3.44) 

9 

3.42 

2-3 

(5.50) 

9-8 

5.14 

3 

(5.75) 

8 

5.89 

3-4 

(5.00) 

8-7 

5.18 

tm 

(8.29) 

8.34 

■ 

4-5 

(4.69) 

7-6 

4.92 

(11.00) 

Mean 
Std.  Dev. 

(5.24) 

(0.48) 

Mean 
Std.  Dev. 

(5.07) 

(0.12) 

Detonation 

Velocity® 

5.08 

±0.14 

a 


b 


Only  gages  6-9  were  used  to  characterize  detonation  process  in  NPN. 
Values  relating  to  other  gages  are  shown  in  parentheses.  See  text  for 
explanation. 

Times  are  for  midpoint  of  leading  edges  of  compressive  waves 
measured  from  a  common  arbitrary  zero  time. 


The  histories  from  both  halves  of  the  target  (all  the  records)  were 
shown  together  in  Figure  23.  The  same  apparently  steady  flows  were 
recorded  in  each  half  at  the  final  two  gage  stations.  In  the  upper 
half,  however,  the  flow  apparently  builds  in  the  region  spanned  by  the 
first  three  gages,  whereas  in  the  lower  half  the  flow  is  steady  at  all 
gage  locations  except  the  first  where  it  is  anomalously  large.  The 
shock  velocity  between  adjacent  gages  is  much  more  uniform  in  the  lower 
half  than  in  the  upper  half,  even  if  gage  1  in  the  upper  half  is 
excluded.  The  shock  velocity  in  the  lower  half  of  the  target  is  consis¬ 
tent  with  steady  flow  at  5.08  km/s. 

Data  Selection 

Histories  1  and  10,  recorded  by  the  longest  and  most  fragile  gages, 
are  anomalous  by  all  three  comparisons.  Because  they  do  not  reproduce 
each  other  and  neither  provides  Information  on  the  steady  detonation 
process  in  NPN,  these  histories  are  not  Included  in  the  Lagrange  analy¬ 
sis  . 

Histories  2-5  from  the  upper  half  of  the  target  show  an  unexpected 
flow  structure  in  which  an  apparent  buildup  occurs  over  the  first  gages 
and  steady  state  is  not  reached  until  the  third  gage  location.  We 
attribute  this  to  a  target  asymmetry  resulting  from  the  filling  proce¬ 
dure  described  previously.  In  the  final  step  of  the  filling  procedure, 
NPN  was  added  through  the  auxiliary  filling  hole  in  the  side  of  the 
chamber.  During  this  step,  voids  developed  in  regions  of  the  target 
chamber  near,  but  not  directly  under,  the  filling  tube.  Attempts  to 
flow  NPN  into  these  regions  were  only  partially  successful,  effectively 
decreasing  the  NPN  specimen  diameter  in  the  upper  half  of  the  target. 

We  consider  it  likely  that  the  reduced  particle  velocity  amplitudes 
recorded  by  the  first  several  gages  in  the  upper  half  of  the  target, 
which  are  those  nearest  the  voids  along  the  side  boundary,  are  a  result 
of  rarefactions  originating  from  these  regions.  The  late  arrival  of  the 
shock  at  gage  1  is  consistent  with  this  Interpretation.  Because  of  this 


target  asymmetry,  we  excluded  the  records  from  the  top  half  of  the 
target,  1-5,  from  the  Lagrange  analysis  and  from  the  characterization  of 
the  detonation  process  In  NPN. 

The  histories  recorded  by  gages  9-6  from  the  bottom  half  of  the 
target  appear  to  be  high  quality  records  of  the  detonation  process  In 
NPN.  These  records  were  therefore  selected  as  the  nonredundant  data  set 
for  flow  characterization  and  Lagrange  analysis.  The  following  discus¬ 
sions  are  based  on  these  records.  As  Indicated  In  Table  3,  the  NPN 
detonation  velocity  calculated  from  particle  velocity  histories  9-6  Is 
5.08  km/s. 

Flow  Characterization 

Histories  6-9  are  shown  superimposed  In  Figure  27,  (The  histories 
were  superimposed  by  adjusting  the  zero  times  for  records  6,  7,  and  8  so 
that  their  "midpoint"  shock  arrival  times  were  the  same  as  that  given 
for  record  9  in  Table  3.)  The  profiles  are  the  same  to  well  within  the 
estimated  experimental  error,  confirming  the  previous  hypothesis  of 
steady  flow.  Although  history  6  shows  a  slight  Increase  above  the 
others  beginning  about  3  ps  after  shock  arrival,  we  do  not  attach 
physical  significance  to  It  because  it  does  not  occur  on  the  redundant 
gage  record  5  nor  on  the  other  histories. 

History  9  lies  within  the  envelope  defined  by  the  superimposed 
profiles  In  Figure  27  and  therefore  Is  a  convenient  representation  of 
the  steady-state  particle  velocity  profile.  The  steady-state  profile 
can  be  idealized  as  consisting  of  three  parts:  a  shock  compression  to  a 
peak  amplitude  of  2.4  kra/s  In  about  0.2  microseconds,  a  relatively  rapid 
decay  to  1.4  kra/s  occurring  In  about  2  ps,  and  then  a  much  slower, 
nearly  linear  decay,  to  1.2  kra/s  in  about  6  ps. 

The  decay  rate  Is  so  slow  In  the  third  region  that  an  alternative 
Idealization  is  to  consider  the  particle  velocity  within  It  to  be  con¬ 
stant  at  1.3  km/s.  The  significance  of  these  alternative  Interpreta¬ 
tions  Is  considered  during  presentation  of  the  results  of  the  thermo¬ 
dynamic  calculations. 


The  experimental  results  may  now  be  summarized.  NPN  In  30-cm- 
dlameter  charges  exhibits  a  steady  detonation  wave  propagating  at  5.08 
±0.1  km/s.  The  Idealized  steady-state  particle  velocity  profile  con¬ 
sists  of  a  shock  compression  to  2.4  km/s  In  about  0.2  ps,  a  decay  to 
about  50%  of  the  peak  amplitude  In  about  2  ps,  and  a  very  slow  further 
decay  or  constant  state  until  the  termination  of  the  experimental  data 
about  6  ps  later.  The  steady  profile  Is  well  represented  by  particle 
velocity  history  9. 

Lagrange  Analysis 

Because  the  flow  Is  steady,  the  Lagrange  analysis  can  be  performed 
by  application  of  the  Ranklne-Hugonlot  shock-jump  equations  expressing 
conservation  of  momentum  and  mass: 

p  =  [Du]/vq  (10) 

V  =  vo[l  -  u/D]  (11) 

where  p  Is  pressure,  D  Is  detonation  velocity,  u  Is  particle  velocity,  v 
Is  specific  volume  (reciprocal  density),  and  the  subscript  0  denotes  the 
Initial  state  Into  which  the  steady  wave  propagates. 

The  Lagrange  pressure  history  In  the  steady-state  region  was  calcu¬ 
lated  by  applying  Equation  (10)  with  D  =  5.08  km/s  and  Vq  =  0.99  m^/Mg 
to  particle  velocity  history  9.  The  result  Is  shown  In  Figure  28.  The 
p-u  and  p-v  paths  are  the  well-known  Rayleigh  lines,  which  may  be  repre¬ 
sented  analytically  from  Equations  (10)  and  (11)  as: 

p/u  =  [D/vg]  (12) 

2 

p/(v  -  vq)  =  (D/vq! 

or,  (13) 

p/(V  -  1)  »  d2/vo 
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Figure  28.  Lagrange  pressure  history  calculated  from  particle  velocity  history, 
gage  9,  second  NPN  experiment. 


where  V  Is  relative  volume,  V  =  v/vq.  in  the  remainder  of  this  section 
it  is  convenient  to  give  volumes  in  terms  of  relative  volume  and  pres¬ 
sures  in  Mhar  (1  Mbar  =  100  GPa),  following  normal  JWL  conventions, 
since  the  results  will  be  used  to  generate  JWL  parameters* 

The  p-V  path  (Rayleigh  line)  calculated  for  PJPN  from  the  second 
experiment  is  shown  in  Figure  29.  This  figure  summarizes  the  experi¬ 
mental  results. 

4.4  TIGER  CALCULATIONS 

Calculations  were  performed  with  the  thermodynamic  equilibrium  code 
TIGER  to  Interpret  and  extend  the  experimental  results  in  Figure  29. 
Specifically,  the  TIGER  calculations  were  used  to  Investigate  the 
ideality  (see  following  paragraphs)  of  the  NPN  reaction  process,  to 
calculate  the  CJ  state,  and  to  extend  the  experimentally  determined 
release  path  to  relative  volumes  approaching  1000.  The  TIGER  calcula¬ 
tions  used  both  the  BKW  and  the  JCZ3  equations  of  state  for  gaseous 
detonation  products.  These  two  equations  of  state  will  be  shown  to 
suggest  rather  different  interpretations  of  the  detonation  process  in 
NPN. 

The  term  ideality  is  commonly  used  to  refer  to  two  properties  of 
detonations;  (1)  charge  size  dependence  and  (2)  completeness  of  the 
chemical  reaction  supporting  steady  detonation.  Consider  first  charge 
size  dependence.  A  real,  quasi-l-D  detonation  process  is  said  to  be 
nonideal  if  the  observed  detonation  parameters  such  as  CJ  pressure  and 
detonation  velocity  are  not  essentially  Independent  of  charge  size. 

Size  dependence  occurs  when  2-D  (lateral)  effects  within  the  reaction 
zone  are  so  significant  that  the  reaction  zone  is  not  well  represented 
as  a  plane  wave.  This  type  of  nonideality  is  believed  to  occur  in  all 
reactive  materials  when  the  ratio  of  charge  size  to  reaction  zone  length 
is  small  enough  and  to  become  negligible  when  the  ratio  is  large 
enough.  In  the  latter  case  we  say  that  "infinite  diameter"  behavior  is 
attained.  The  present  TIGER  calculations  are  Inherently  1-D  and  can  be 
sensibly  compared  only  with  results  from  "infinite  diameter"  experi¬ 
ments  . 
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Figure  29.  NPN  pressure-relative  volume  compression  and  release  paths, 
second  NPN  experiment. 

The  locus  of  Lagrange  p-V  states  calculated  from  particle 
velocity  records  6-9. 


Some  reactive  materials  such  as  NPN  have  longer  reaction  zones  than 
conventional  military  explosives.  These  materials  are  sometimes  called 
"nonldeal"  because  they  do  not  exhibit  Infinite  diameter  behavior  except 
for  relatively  large  charges.  In  fact,  like  other  reactive  materials, 
they  may  support  either  Ideal  or  nonlieal  detonations  depending  on 
charge  size  and  confinement.  In  the  present  work,  we  assume  that  the 
0.3-m-dlaraeter  contained  NPN  charges  studied  here  are  large  enough  to 
exhibit  "Infinite  diameter"  behavior.  This  assumption  Is  reasonable 
because  ANFO  data  reported  by  Finger  et  al.^®  Indicates  that  this  would 
be  a  good  approximation  for  ANFO  and  because  NPN  Is  expected  to  have  a 
shorter  reaction  zone  because  of  Its  finer  particle  size,  greater 
Intimacy  of  mixing,  and  enhanced  detonation  pressures  and  tempera- 
tures.  The  NPN  results  reported  by  Franzen  and  Wlsotskl^^  are  also 
consistent  with  this  assumption.  We  therefore  Interpret  the  results  of 
the  present  NPN  experiments  as  Infinite  diameter  behavior  and  do  not 
consider  this  type  of  nonideality  further. 

The  other  detonation  property  referred  to  by  the  term  Ideality  Is 
the  completeness  of  the  chemical  reaction  supporting  the  detonation. 
Kinetic  constraints  may  prevent  the  completion  of  the  chemical  reaction 
minimizing  free  energy  and  the  attainment  of  thermodynamic  equilibrium 
at  the  back  of  the  steady-state  reaction  zone.  Then,  as  In  the  case  of 
size  dependent  reactions,  not  all  the  thermodynamically  predicted  chemi¬ 
cal  energy  Is  available  to  support  the  1-D  detonation  and  the  process  Is 
called  nonldeal.  Certain  types  of  such  "chemical  nonideality"  can  be 
accounted  for  In  thermodynamic  calculations  If  the  specific  reaction 
phenomenology  resulting  from  the  constraints  Is  known. 

In  the  remainder  of  this  report  we  use  the  term  Ideality  to  refer 
to  the  effects  of  chemical  reaction  constraints  rather  than  to  size 
dependence.  By  this  definition,  a  detonation  process  Is  Ideal  If  the 
experimentally  measured  and  the  thermodymanlcally  calculated  values  for 
the  CJ  state  agree  within  the  accuracies  of  the  measurement  and  the 
calculation. 


As  the  first  step  In  Interpreting  the  experimental  results  using 
BKW,  we  calculated  the  Ideal  NPN  CJ  state  and  CJ  release  Isentrope  using 
the  SRI  BKW  parameters  for  the  gaseous  species  and  rejecting  liquid 
water.  The  resulting  CJ  state  values  were  D  =  6.33  km/s,  p  =  10.9  GPa, 

V  =  0.724  m^/Mg,  and  u  =  1.70  km/s.  Figure  30  shows  the  p-V  repre¬ 
sentations  of  this  CJ  state  and  the  corresponding  CJ  release  Isentrope 
to  a  relative  volume  of  about  2  (the  complete  calculated  CJ  Isentrope  Is 
shown  later  in  a  log-log  plot).  Figure  30  Includes  the  experimentally 
determined  compression  and  release  states  for  comparison  with  the  Ideal 
BKW  calculations. 

The  BKW  value  for  D  Is  significantly  greater  than  the  experimen¬ 
tally  observed  steady  propagation  velocity,  so  the  calculated  CJ  state 
does  not  lie  on  the  experimental  Rayleigh  In  Figure  30.  Thus,  on  the 
bases  of  BKW/TIGER  calculations  and  our  definition  of  nonideality,  the 
detonation  process  In  NPN  must  be  nonideal. 

Assuming  that  the  detonation  process  In  NPN  is  nonideal  and  that 
the  Lagrange  gage  experiment  does  not  define  the  end  of  the  Initial 
steady-state  reaction  zone,  we  must  postulate  a  specific  calculable  non¬ 
ideal  reaction  phenomenology  to  determine  a  final  detonation  state  and 
subsequent  adiabatic  expansion  states  in  NPN  with  BKW/TIGER.  The 
present  experiments  and  the  other  available  NPN  Information  provide  no 
guidance  because  they  do  not  monitor  the  detonation  product  behavior  In 
the  regime  where  the  nonideal  behavior  occurs. 

'll  OO 

For  ANFO,  however,  Johnson  et  al.-"^  and  Mader  have  reported  the 
results  of  underwater  cylinder  expansion  experiments  (aquarium  tests) 
and  have  presented  methods,  based  on  a  specific  nonideal  reaction  pheno¬ 
menology,  for  numerically  modeling  such  expansions  with  thermodynamic 
BKW  calculations.  In  the  present  work,  we  assuraraed  that  the  same  quali¬ 
tative  phenomenology  also  applies  to  NPN.  Our  justification  for  this 
assumption  Is  that  we  expect  ANFO  and  NPN  to  react  similarly  and  that 
the  aquarium  teat  provides  ANFO  product  expansion  data  to  volumes  large 
enough  to  provide  a  relevant  check  against  the  predictions  based  on  the 
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postulated  phenomenology  and  BKW  calculations.  Following  Johnson  et  al. 
and  Mader,  we  made  the  following  assumptions  and  calculations  to  model 
the  nonideal  detonation  In  NPN. 


First,  we  assumed  that  only  a  portion  of  the  ammonium  nitrate 
reacts  to  support  the  steady-state  detonation  and  determined  the  amount 
of  AN  that  must  be  held  Inert  to  match  the  steady  detonation  velocity 
observed  In  the  Lagrange  gage  experiments.  We  calculated  a  CJ  detona¬ 
tion  state  based  on  this  fraction  of  Inert  AN.  Next  we  assumed  that.  In 
the  Initial  part  of  a  real  adiabatic  expansion  from  this  state  (such  as 
might  be  observed  In  a  very  large  Lagrange  gage  experiment),  the  remain¬ 
ing  AN  would  be  consumed  irreversibly.  Johnson  et  al.  and  Mader  report¬ 
ed  that  the  numerical  simulations  of  the  ANFO  aquarium  tests  were  insen¬ 
sitive  to  the  assumptions  used  In  modeling  the  secondary  AN  reaction. 

We  postulated  a  specific  Irreversible  process  and  calculated  the  cor¬ 
responding  p-V  path  and  entropy  Increase.  Finally,  we  assumed  that  the 
remainder  of  the  expansion  to  atmospheric  pressure  that  would  be 
observed  In  an  appropriate  experiment  would  occur  Isentroplcally,  and  we 
calculated  the  corresponding  Isentrope.  Details  of  the  NPN  calculations 
are  given  below. 

We  first  calculated  nonideal  CJ  detonation  velocities  for  various 
fractions  of  Inert  AN.  The  Inert  AN  was  modeled  as  simply  as  possible 
within  the  options  available  In  the  TIGER  solid  equation  of  state.  The 
AN  was  assumed  to  be  Incompressible  and  to  have  a  constant  specific  heat 
=  32.82  cal/degree  mole,  a  standard  state  heat  of  formation  =  -87.3 
kcal/mole,  and  a  standard  state  entropy  =  33.8  cal/degree  mole  (esti¬ 
mated  from  values  for  related  compounds).  In  these  and  the  subsequent 
TIGER  calculations  for  NPN,  the  methocel,  which  Is  only  0.2  wt%  of  the 
NPN  and  is  a  complex  C-H-0  molecule,  was  simply  accounted  for  by 
increasing  the  amount  of  methyl  alcohol  by  0.2  wt%. 

The  experimental  detonation  velocity  of  5.08  km/s  was  matched  by 
BKW  calculations  with  51.5  wt%  AN  Inert.  Thus  a  nonideal  CJ  detonation 
with  51.5  wt%  AN  inert  reproduces  the  experimentally  determined  Rayleigh 
line.  The  corresponding  CJ  state  variables  determined  from  the  TIGER 
calculations  are  p  =  5.44  GPa,  u  »  1.08  km/s,  and  v  »  0.780  Mg/m^. 
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The  nonideal  BKW/TIGER  CJ  pressure  and  particle  velocity  are  less 
than  the  final  values  from  the  second  NPN  experiment,  as  expected.  This 
Is  consistent  with  the  Interpretation  that  the  flow  observed  In  this 
experiment  Is  entirely  within  the  primary  AN  reaction  zone.  The  calcu¬ 
lated  nonideal  CJ  detonation  state  then  defines  the  end  of  the  reaction 
zone. 

We  can  now  estimate  a  reaction  time  for  the  primary  AN  reaction  by 
assuming  the  first  steady-state  particle  velocity  profile  Idealization 
discussed  In  the  previous  section.  In  this  Idealization,  the  final 
portion  of  the  history  Is  taken  as  a  slow  linear  decay.  We  extrapolated 
the  approximate  decay  rate  to  the  calculated  CJ  state  to  estimate  a 
primary  AN  reaction  time  of  11  ±3  ps. 

We  noted  In  Section  4.2  that  the  last  part  of  the  experimental 
particle  velocity  history  may  also  be  Idealized  as  a  constant  state. 

With  this  Interpretation,  the  reaction  time  estimate  above  Is  clearly 
Invalid.  However,  the  constant  state  Interpretation  Is  not  believed  to 
be  consistent  with  the  present  BKW  calculations  for  the  following 
reason. 

The  only  plausible  explanation  we  can  provide  for  a  steady  wave 
followed  by  a  truly  constant  state  Is  that  the  NPN  is  overdriven.  (This 
would  occur  If  a  constant  particle  velocity  greater  than  the  CJ  particle 
velocity  were  applied  to  the  rear  boundary  of  the  NPN;  the  driver  Input 
criteria  discussed  In  the  Section  4.2  suggest  that  this  Is  extremely 
unlikely.)  However,  if  NPN  Is  overdriven,  then  Its  non-overdrlven 
steady-state  detonation  velocity  is  less  than  our  experimental  value  of 
5.08  kra/s.  But  for  the  BKW/TIGER  calculations  to  be  consistent  with 
such  a  detonation  velocity,  essentially  all  the  AN  would  have  to  be  held 
Inert.  Because  we  consider  this  to  be  an  unreasonable  condition,  on  the 
basis  of  these  BKW  calculations  we  conclude  that  the  NPN  Is  apparently 
not  overdriven,  and  In  this  section  we  discard  the  possibility  of  Ideal¬ 
izing  the  last  part  of  the  experimental  record  as  a  constant  state. 

To  calculate  the  expansion  from  the  nonideal  CJ  state  to  the  state 
where  the  remaining  AN  Is  consumed,  we  must  define  a  specific  secondary 
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reaction  process.  From  the  results  of  Johnson  et  the  calcula- 

tlonal  simulations  of  the  aquarium  tests  are  not  sensitive  to  the  speci¬ 
fic  process  chosen.  We  chose  to  model  the  process  as  a  CJ  deflagration 
centered  on  the  nonideal  detonation  state. 

This  choice  apparently  Implies  a  second  steady  reactive  flow 
traveling  slower  than  the  steady  detonation  observed  experimentally, 
with  a  constant  state  of  growing  duration  between  them.  A  constant 
state  Is  compatible  with  our  flow  idealization.  However,  a  growing 
constant  state  Is  not  physically  realistic  because  the  later  a  particle 
enters  the  flow,  the  longer  it  will  wait  to  undergo  the  secondary  reac¬ 
tion  even  though  it  Is  preconditioned  identically  to  particles  entering 
the  flow  earlier.  Nevertheless,  the  BKW/TIGER  release  path  is  Intended 
to  be  used  only  as  an  intermediate  step  In  generating  a  JWL  isentrope 
through  the  CJ  state.  Because  the  JWL  Isentrope  form  always  has  a 
positive  second  derivative  with  respect  to  relative  volume,  the  physi¬ 
cally  unrealistic  hydrodynamic  behavior  will  not  be  present  in  a  JWL  EOS 
fitted  to  the  nonldeal  BKW  CJ  isentrope. 

The  p-V  path  for  the  postulated  deflagration  process  is  the 
Rayleigh  line  between  the  nonldeal  detonation  state  and  the  CJ  deflagra¬ 
tion  state.  From  BKW/TIGER  calculations,  at  the  CJ  deflagration  state: 
p  =  2.96  GPa,  V  *  1.19  m^/Mg,  and  S  =  1.727  cal/Kg,  =  1.292  cal/Kg), 
where  S  is  specific  entropy  and  the  subscript  1  refers  to  the  nonldeal 
detonation  state.  The  remainder  of  the  release  path  is  modeled  simply 
as  the  BKW  Isentrope  through  this  CJ  deflagration  state. 

The  nonldeal  BKW  adiabatic  release  path  generated  as  described 
above  to  model  and  extend  the  results  of  the  second  NPN  experiment  then 
consists  of  the  Rayleigh  line  from  the  nonldeal  CJ  detonation  state  for 
51.5  wt%  AN  inert  to  an  appropriately  centered  ideal  (AN  100%  active)  CJ 
deflagration  state  and  an  ideal  BKW  Isentrope  through  this  deflagration 
state.  Figure  31  shows  the  nonldeal  BKW  detonation  state  and  the  adlab- 
tlc  release  path  to  a  relative  volume  of  about  2.  Figure  31  also 
includes  the  experimental  data  and  the  Ideal  BKW  CJ  isentrope  for  com¬ 
parison. 
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Figure  32  shows  the  CJ  states  and  full  release  paths  to  the  largest 
relative  volumes  calculated  for  the  ideal  and  nonideal  cases  on  a  log- 
log  scale.  (The  isentrope  calculations  were  carried  out  to  the  largest 
volume  expansions  possible;  the  calculations  terminated  when  tempera¬ 
tures  dropped  below  the  limits  built  into  the  EOS.)  The  ideal  BKW  CJ 
isentrope  and  the  lower  portion  of  the  nonideal  BKW  CJ  Isentrope  are 
indistinguishable  in  this  figure  because  the  entropy  productions  in  the 
one-step  and  two-step  AN  reaction  processes  are  nearly  the  same.  The 
significance  of  these  thermodynamic  calculations  is  discussed  after  the 
JCZ3  calculations  are  presented. 

JCZ3  Calculations 

As  the  first  step  in  interpreting  the  experimental  results  using 
JCZ3,  we  again  calculated  the  ideal  NPN  CJ  state  and  CJ  release  Isen¬ 
trope.  We  used  the  SRI  JCZ3  parameters  for  the  gaseous  species,  reject¬ 
ed  liquid  water,  and  assumed  ideal  detonation.  The  resulting  CJ  state 
values  were  D  =  5.48  km/s,  p  =  7.93  GPa,  v  *  0.731  ra^/Mg,  and  u  =  1.43 
kra/s.  This  CJ  state  and  the  ideal  CJ  release  isentrope  to  a  relative 
volume  of  about  1.8  are  shown  in  Figure  33  along  with  the  experimental 
Rayleigh  line  presented  previously.  (The  full  release  Isentrope  is 
shown  later.) 

In  sharp  contrast  to  the  BKW  calculations,  the  JCZ3  CJ  detonation 
calculation  with  AN  fully  active  reproduces  the  experimentally  observed 
detonation  velocity  within  the  estimated  accuracies  of  the  calculation 
and  experiment.  (The  experimental  measurement  was  previously  estimated 
to  have  an  uncertainty  of  iO.l  km/s;  the  calculations  are  judged  to  have 
an  uncertainty  of  about  ±5%.)  Therefore,  by  the  definition  given 
previously,  with  respect  to  the  JCZ3  equation  of  state,  the  detonation 
process  in  NPN  is  ideal. 

The  conclusion  that  the  detonation  process  observed  in  NPN  in  the 
Lagrange  gage  experiment  is  ideal  with  respect  to  JCZ3  follows  from  our 
definitions  and  is  convenient  for  indicating  that  the  ideal  JCZ3  cal¬ 
culations  much  more  closely  replicate  the  observed  NPN  behavior  than  the 
ideal  BKW  calculations.  It  could,  of  course,  be  argued  that  the 
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Figure  33.  NPN  JCZ3  pressure-relative  volume  CJ  state  and  partial  release 
isentrope,  NPN  ideal. 

The  experimental  results  (solid  line)  from  Figure  29  are  included 
for  comparison. 
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discrepancy  is  significant  and  the  NPN  detonation  process  is  nonideal. 

In  fact,  in  the  remainder  of  this  section,  to  proceed  toward  our  primary 
objective  of  developing  initial  JWL  parameters  for  NPN,  we  Improve  the 
correlation  between  experiment  and  calculation  by  treating  the  NPN 
detonation  process  as  nonideal. 

To  obtain  precise  agreement  in  CJ  detonation  velocity  between  the 
JCZ3  calculations  and  the  experiment,  we  again  assumed  that  not  all  the 
AN  reacted  to  support  the  steady  flow.  Ue  determined  that  agreement  was 
obtained  if  14.5  wt%  of  the  AN  was  held  inert  in  the  calculations.  The 
corresponding  CJ  state  variables  are  p  =  6.48  GPa,  u  *  1.26  km/s,  and 
V  =  0.744  Mg/m^.  We  take  this  to  be  the  nonideal  JCZ3  CJ  state  in  NPN. 


The  calculated  nonideal  JCZ3  CJ  state  is  reached  during  the 
Lagrange  gage  NPN  experiment,  in  distinction  to  the  situation  with  the 
BKU  calculations.  The  JCZ3  values  for  the  nonideal  CJ  pressure  and 
particle  velocity  are  located  on  the  corresponding  experimentally  deter¬ 
mined  histories  either  at  the  end  of  the  first  rapid  rarefaction  or 
early  in  the  linear  decay  region.  On  this  basis,  we  estimate  the  reac¬ 
tion  time  as  2  (+5  -1)  ps  for  the  primary  AN  reaction  (the  reaction  that 
supports  the  observed  detonation  velocity) . 

The  experimentally  observed  steady  flow  persists  to  lower  pressures 
and  particle  velocities  than  the  primary  AN  reaction,  which  terminates 
at  the  nonideal  JCZ3  CJ  point.  This  suggests  either  that  (1)  some 
secondary  reaction  occurs  near  the  CJ  state,  maintaining  the  steady  flow 
and  causing  the  slow  pressure  and  particle  velocity  decay  rates,  or  that 
(2)  the  flow  is  not  rigorously  steady  below  the  nonideal  CJ  state, 
although  we  failed  to  resolve  the  dispersion  from  our  measurements.  In 
the  latter  case  we  would  still  have  to  account  for  the  unusually  low 
decay  rates  of  the  Lagrange  stress  and  particle  velocity  histories. 

This  slow  decay  differs  from  the  behavior  we  have  observed  In  other 
explosives  and  is  an  interesting  result.  We  consider  these  points 
further  to  calculate  the  JCZ3  nonideal  adiabatic  CJ  release  path. 


In  calculating  a  JCZ3  nonideal  CJ  release  path,  we  must  again 
postulate  a  calculable  reaction  phenomenology  as  we  did  in  the  case  of 
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the  BKW  calculations.  We  choose  not  to  use  the  ANFO  BKW  phenomenology 
Itself,  however,  because  It  was  developed  and  validated  specifically  for 
BKW  calculations.  To  apply  the  method  to  JCZ3  calculations,  we  would 
have  had  to  arbitrarily  choose  a  secondary  reaction  process  and  a  final 
pressure  for  terminating  the  process  with  no  means  of  Justifying  the 
choice. 

We  calculated  JCZ3  release  Isentropes  for  two  assumed  processes. 

In  the  first,  the  AN  that  did  not  support  the  steady  detonation  was 
assumed  also  to  remain  Inert  In  the  subsequent  adiabatic  release.  This 
assumption  Is  computationally  convenient  and,  because  most  (85.5  wt%)  of 
the  AN  Is  consumed  In  the  primary  reaction,  may  be  sufficiently 
accurate.  It  falls  to  account,  however,  for  the  anomalously  slow  decay 
rates  observed  In  the  experimental  records. 

The  second  assumed  release  process  was  more  complex  and  was 
motivated  by  the  shape  of  the  experimentally  observed  pressure  and 
particle  velocity  histories.  In  developing  our  second  postulated  pro¬ 
cess  from  the  nonideal  CJ  state,  we  were  motivated  by  the  experimental 
observation  that  the  pressure  decay  rate  was  so  slow  that  the  pressure 
could  be  Idealized  as  being  constant.  We  now  attribute  this  not  to  an 
overdriven  detonation  (which  we  consider  highly  unlikely  because  of  the 
Initiator  selection  criteria  described  earlier),  but  to  a  secondary  AN 
reaction.  Furthermore,  If  the  particle  velocity  Is  also  Idealized  as 
constant,  then  the  particle  velocity  gradient  Is  zero.  Therefore,  by 
the  continuity  condition  for  uniaxial  flow  (Equation  3,  Section  2),  the 
specific  volume  Is  constant.  Thus,  the  experimental  profiles  suggest  a 
secondary  reaction  occurring  at  nearly  constant  p  and  V.  We  therefore 
postulated  that  (1)  such  a  reaction  exists,  (2)  It  consumes  the  remain¬ 
ing  AN,  and  (3)  the  constant  p  and  V  values  are  those  calculated  for  the 
nonideal  CJ  state.  We  do  not  argue  the  plausibility,  or  even 
possibility,  of  such  a  process,  but  only  state  that  It  Is  consistent 
with  the  experimental  results. 

We  calculated  JCZ3  Isentropes  for  both  postulated  processes.  In 
the  first  case,  the  JCZ3  Isentrope  Is  the  frozen  Isentrope  for  1A.5  wt% 


AN  Inert  through  the  nonideal  CJ  state.  In  the  second  case.  It  Is  the 
Ideal  (all  AN  active)  Isentrope  through  the  nonideal  CJ  state.  These 
two  Isentropes  plus  the  Ideal  JCZ3  CJ  isentrope  are  shown  In  their 
entirety  In  Figure  34.  Note  that  the  JCZ3  calculations  extend  to  rela¬ 
tive  volumes  approaching  1000. 

We  selected  the  nonideal  CJ  release  isentrope  based  on  the  second 
postulated  process  as  the  preferred  JCZ3  representation  for  the  NPN 
steady  detonation  and  expansion.  This  path,  to  relative  volumes  of 
about  1.8,  is  shown  in  Figure  35  along  with  the  Rayleigh  line  from  the 
experiment  and  the  ideal  JCZ3  CJ  state  and  release  isentrope.  We  used 
this  path  to  develop  initial  JWL  parameters  for  NPN. 

4.5  CALCULATION  OF  JWL  COEFFICIENTS 

Initial  JWL  coefficients  were  developed  for  NPN  using  the  proce¬ 
dures  described  in  Section  2.  The  nonideal  JCZ3  CJ  state  and  release 
path  (Figure  34,  dashed  curve,  and  Figure  35)  were  selected  as  the  basis 
characterization  data  for  use  in  developing  the  JWL  parameters.  We  had 
intended  to  also  use  the  BKW  nonideal  CJ  adiabatic  release  path  for  this 
purpose.  However,  the  failure  of  the  calculations  to  reach  volume 
expansions  large  enough  to  define  realistic  low  pressure  limits  for  the 
slope  of  the  CJ  release  adlabat  prevented  this. 

The  fitting  procedures  produced  JWL  parameters  consistent  with  the 
nonideal  JCZ3  CJ  release  adlabat  as  well  as  with  the  three  constraints 
on  CJ  pressure,  slope,  and  energy.  Equations  (7),  (8),  and  (9).  The 
only  choice  required  was  a  value  of  Vj..  We  selected  =  25  and  encoun¬ 
tered  no  difficulties  in  completing  the  fitting  process. 

The  JWL  parameters  calculated  for  NPN  are  listed  in  Table  4  in 
conventional  JWL  format  and  units.  Note  that  in  this  format  the  CJ 
relative  volume  is  determined  from  F.  For  convenience  the  values  of  the 
various  CJ  state  flow  parameters  not  given  in  Table  4  are  summarized 
here:  u  *  0.126  cm/ps,  v  =  0.744  m^/Mg,  and  V  =  0.7514. 
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Figure  34.  Three  JCZ3  calculations  of  NPN  pressure-relative  volume  CJ  state 
and  release  isentrope. 
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Figure  35.  NPN  JCZ3  pressure-relative  volume  CJ  state  and  partial  release 
adiabat,  NPN  nonideal. 

The  experimental  results  (light  solid  line)  from  Figure  29  and  the 
ideal  JCZ3  results  (dashed  line)  from  Figure  33  are  included  for 
comparison.  See  text  for  discussion  of  nonideal  reaction  process. 
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Table  4.  NPN  JWL  coefficients. 


Inlclal 

State 

CJ  State 

Equatlon-of-State  Parametera 

Po 

(Mg/m^) 

P 

(Mbar) 

D 

(cm/)is) 

Eo 

(Mbar) 

■ 

A 

(Mbar) 

B 

(Mbar) 

C 

(Mbar) 

B 

R2 

U) 

1.01 

0.0648 

0.508 

0.3187 

3.023 

6.783 

0.1824 

0.007212 

7.628 

2.305 

0.345 

The  fitted  (JWL)  and  basis  (JCZ3)  CJ  isentropes  are  shown  in  the 
p-V  plane  in  Figure  36.  The  JWL  and  JCZ3  adiabatic  expansion  coeffi¬ 
cients  Y»  where  y  “  (d  In  p/d  In  V)g,  for  these  paths  are  shown  in 
the  plane  in  Figure  37.  No  adjustments  are  made  in  our  fitting 
method  to  affect  the  agreement  between  the  curves  in  Figure  37. 

4.6  DISCUSSION 

The  NPN  JWL  parameters  in  Table  4  are  the  major  result  of  the 
effort  described  in  this  section.  They  incorporate  both  experimental 
observations  of  compression  and  release  behavior  in  the  reaction  zone  as 
well  as  compatible  thermodynamic  calculations  based  on  JCZ3,  a  state-of- 
the-art  gaseous  products  equation  of  state.  These  .JWL  parameters  are 
judged  appropriate  for  use  in  Initial  hydrocode  calculations  to  evaluate 
NPN  behavior  in  simulators  for  DNA  applications. 

Most  of  the  specific  results  of  the  NPN  characterization  subtasks 
were  discussed  as  the  results  were  presented.  Here  we  briefly  consider 
some  general  implications  of  this  work. 

The  pretest  effort  to  develop  procedures  for  obtaining  uniform  NPN 
at  densities  around  1  Mg/ra^  Indicated  the  difficulty  in  keeping  NPN 
uniform  at  these  low  densities.  We  conclude  that  delicate  NPN  handling 
procedures  would  be  necessary  to  get  uniform  performance  in  simulators 
from  low  density  NPN.  However,  in  DNA  applications,  as  in  the  mining 
applications  for  which  NPN  was  developed.  It  is  likely  that  NPN  would  he 
used  at  higher  densities  to  enhance  output.  In  such  cases,  tamping  the 
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Figure  37.  Gamma  versus  relative  volume  for  NPN  JWL  CJ  release  isentrope 
and  NPN  CJ  release  path  being  fitted. 

The  two  gammas  are  not  equal  at  the  CJ  relative  volume  because 
the  NPN  CJ  release  path  is  not  tangent  to  the  Rayleigh  line  at 
this  point. 
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charge  during  loading  to  achieve  maximum  Irreversible  compaction  every¬ 
where  may  result  In  a  uniform  as  well  as  a  high  output  material.  The 
tendency  of  NPN  to  compact  Irreversibly  should  be  borne  In  mind  when 
considering  It  for  a  specific  application. 

Our  experiments  showed  that  NPN  supported  a  steady  detonation  with 
a  propagation  velocity  of  5.08  km/s  In  confined  0.3-m-dlaraeter 
charges.  From  ANFO  results,  we  Inferred  that  this  value  Is  a  good 
approximation  to  the  Infinite  diameter  velocity.  This  Inference  Is  also 
compatible  with  the  experimental  results  of  Franzen  and  Wlsotskl,  for  a 
lOO-pound  NPN  charge^®  If  the  propagation  velocities  they  report 
(between  2.8  and  4.1  km/s)  are  Interpreted  as  being  from  the  buildup 
region.  There  Is  an  Implication  In  the  Introduction  to  Reference  20 
that  these  values  represent  steady  detonation;  however,  from  reviewing 
the  data  and  consulting  with  the  authors,  we  determined  that  this  Impli¬ 
cation  was  unintentional. 

The  Inference  that  our  Lagrange  gage  experiment  approximates 
Infinite  diameter  behavior  Is  also  qualitatively  compatible  with 
Mallory's  greater  value  of  5.206  km/s  for  NPN  of  higher  density  (1.298 
Mg/m^),  In  1-lnch-dlameter  heavily  confined  columns.  However,  because 
the  difference  In  detonation  velocities  Is  rather  small  considering  the 
significant  density  difference  between  Mallory's  material  and  ours.  It 
Is  possible  that  one  or  both  of  the  velocities  do  not  represent  Infinite 
diameter  behavior.  Nevertheless,  In  view  of  the  large  diameter  of  our 
charges  and  the  supporting  ANFO  Information,  we  believe  that  5.08  km/s 
Is  a  good  approximation  to  the  Infinite  diameter  detonation  velocity  In 
NPN  of  1.01  Mg/m^  Initial  density. 

The  Lagrange  particle  velocity  histories  Indicated  two  unusual 
hydrodynamic  properties  of  the  detonation  process  In  NPN.  First,  the 
flow  was  steady  throughout  the  estimated  1-D  test  time  at  each  gage  (at 
least  5-10  ps  after  shock  arrival).  Second,  the  decay  rates  of  the 
Lagrange  particle  velocity  and  Inferred  pressure  histories  became  very 
slow  or  zero  about  two  microseconds  after  shock  arrival.  We  have  not 
observed  either  of  these  detonation  properties  In  other  HEs  we  have 


105 


studied  with  this  technique.  Assuming  the  results  to  be  reproducible 
and  not  an  artifact  of  the  technique,  we  conclude  that  the  reaction 
phenomenology  In  NPN  Is  significantly  different  from  that  In  the  other 
HEs  we  have  studied. 

We  found  significant  differences  In  the  Implications  of  thermodyna¬ 
mic  TIGER  calculations  modeling  the  experiments  with  the  BKI'f  and  with 
the  JCZ3  equations  of  state.  The  TIGER/BKW  calculations  of  detonation 
velocity  suggested  that  the  detonation  Is  definitely  nonideal,  but  the 
TIGER/JCZ3  calculations  Implied  that  It  Is  only  borderline  nonideal. 

The  TIGER/JCZ3  calculations  Indicated  that  most  (85.5  wt%)  of  the  AN 
reacts  to  support  the  detonation,  whereas  the  TIGER/BKW  calculations 
Implied  that  only  half  the  AN  reacted  to  support  the  detonation.  The 
estimated  nonideal  CJ  state  occurred  during  the  gage  recording  times,  at 
a  pressure  of  6.48  GPa,  according  to  the  JCZ3  calculations,  but  after 
the  gage  recording  times,  at  5.44  GPa,  In  the  BKW  calculations. 

Because  of  the  differences  In  the  BKW  and  JCZ3  nonideal  CJ  calcu¬ 
lations,  we  postulated  different  secondary  reaction  processes  In  the  two 
cases:  a  nearly  constant  volume  reaction  occurring  In  the  region  of  the 
CJ  pressure  for  the  JCZ3  calculations  and  a  steady  deflagration  spread 
over  a  wide  range  of  pressures  and  volumes  for  the  BKW  calculations. 

Both  are  artifacts  to  permit  the  completion  of  the  calculations  and 
provide  a  reasonable  estimate  for  the  locus  of  p-V  states  along  the  NPN 
CJ  adlabat  for  constraining  the  JWL  fit.  Since  each  is  known  to  have 
certain  unrealistic  features,  neither  Is  represented  to  be  the  actual 
process  occurring  In  NPN;  unfortunately,  the  present  work  cannot  indi¬ 
cate  which  Is  preferable.  We  conclude  that  NPN  (and  probably  Its  close 
relative  ANFO)  has  a  complex  unknown  AN  decomposition  phenomenology  that 
should  be  studied  further. 

The  NPN  JWL  parameters  developed  here  were  based  on  the  JCZ3  calcu¬ 
lations  because  the  BKW  calculations  terminated  before  providing  the  low 
pressure  limiting  data  required  for  our  JWL  fitting  method.  We  note 
that  most  of  the  present  JWL  parameters  are  within  the  normal  ranges  for 
such  coefficients  developed  from  cylinder  tests,  with  the  possible 
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exception  of  the  nonlinear  coefficients  and  R2.  The  present  values 
for  these  coefficients  have  a  ratio  typical  of  cylinder  test  coeffi¬ 
cients,  but  are  somewhat  larger  than  published  coefficients.  Larger 
nonlinear  coefficients  have  the  effect  of  bringing  on  the  constant  gamma 
(adiabatic  exponent)  behavior  along  the  CJ  isentrope  at  somewhat  smaller 

relative  volumes.  For  comparison,  we  list  the  cylinder  test  ANFO  JWL 
30 

parameters  and  our  values  In  Table  5  and  plot  the  ANFO  and  NPN  JWL  CJ 
Isentropes  in  Figure  38. 


Table  5.  Comparison  of  NPN  and  ANFO  JWL  coefficients. 


Initial 

State 

CJ  State 

Equatlon-of-State  Parameters 

PQ 

(Mg/m^) 

P 

(Mbar) 

D 

(cm/ns) 

Eo 

(Mbar) 

■ 

A 

(Mbar) 

B 

(Mbar) 

C 

(Mbar) 

B 

B 

0) 

1.01 

0.0648 

0.508 

0.03187 

3.023 

6.783 

0.1824 

0.007212 

7.628 

2.305 

0.345 

ANFO® 

0.850 

- 1 

0.060 

0.465 

0.0325 

2.063 

0.4760 

0.00524 

0.00720 

3.5 

1  0.9 

0.31 

We  conclude  this  subsection  with  a  digression  concerning  the  avail- 

OC 

able  detonation  energy  of  NPN.  The  report  by  McMullan  contains  a 
useful  compilation  of  properties  for  several  HEs  of  interest  to  DNA 
including  NPN.  However,  it  cites  a  value  for  the  available  energy  of 
NPN,  namely,  1180  cal/g,  that  is  unusually  high  compared  with  the  values 
for  similar  HEs.  The  available  energy  value  widely  quoted  for  ANFO,  for 
example,  and  given  in  reference  25,  is  920  cal/g.  We  believe  that  the 
NPN  value  in  Reference  25  is  calculated  on  a  different  basis  than  the 
values  for  the  other  HEs  and  could  therefore  lead  to  erroneous  conclu¬ 
sions  regarding  the  performance  of  NPN  relative  to  these  materials  in 
DNA  applications. 

There  is  no  universally  accepted  definition  for  the  term  available 
energy.  In  practice,  computed  values  for  available  chemical  energy  of  a 
detonating  HE  are  often  either  (1)  maximum  heats  of  detonation 

^From  Reference  30 
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Figure  38.  NPN  (solid)  and  ANFO  JWL  CJ  isentropes. 
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calculated  by  assuming  that  the  initial  and  final  states  are  298  K  and 
101  kPa  and  that  the  detonation  products  are  the  most  stable  ordered  set 
H2O,  CO2,  C(s),  and  N2,  or  (2)  smaller  heats  of  detonation  calculated 
assuming  the  same  initial  and  final  states  but  a  more  realistic  detona¬ 
tion  product  composition  estimated  from  a  thermodynamic  equilibrium  code 
such  as  TIGER,  usually  for  a  relative  expansion  volume  around  10  where 
product  composition  freezeout  is  estimated  to  occur.  In  either  case,  it 
is  necessary  to  specify  whether  the  H2O  is  in  liquid  or  gaseous  form, 
with  the  latter  assumption  being  usual. 

We  performed  both  types  of  calculations  for  NPN.  For  the  baseline 
NPN  composition,  we  calculate  a  maximum  heat  of  detonation  of  907  cal/g 
with  respect  to  H20(g).  This  is  the  value  that  we  believe  should  be 
compared  with  the  other  available  energy  values  in  Reference  25.  (For 
our  composition,  and  using  the  second  definition,  we  calculated  a  heat 
of  detonation  of  877  cal/g  using  the  product  concentrations  on  the  non- 
ideal  TIGER/ JCZ3  CJ  adlabat  at  a  relative  volume  of  about  9.  This  value 
is  larger  than  the  NPN  Eq  value  from  Tables  4  and  5.) 

We  consider  it  likely  that  the  high  available  energy  value  for  NPN 
in  reference  25  is  calculated  with  respect  to  H20(Jl),  whereas  the  values 
for  other  HEs  in  reference  25  are  with  respect  to  H20(g).  When  the 
values  are  calculted  on  a  consistent  basis,  NPN  and  ANFO  appear  to  be 
very  similar  in  available  energy,  with  ANFO  being  slightly  more  ener¬ 
getic  per  gram  of  HE. 
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SECTION  5 


RECOMMENDATIONS 


Results  of  the  present  NA/NP  and  NPN  characterization  efforts  lead 
us  to  make  five  recommendations  to  DNA  In  the  areas  of  HE  use  and 
characterization.  Recommendations  1,  3,  and  S  are  recommended  strongly. 

Recommendation  1:  We  strongly  recommend  that  baseline  NA/NP  be 
considered  for  use  In  applications  calling  for  low  cost  liquid  explo¬ 
sives  similar  In  performance  to  nltromethane. 

NA/NP  is  a  very  promising  family  of  liquid  explosives.  It  Is  a 
uniform,  high  quality  HE  with  performance  very  similar  to  nltromethane, 
but  the  material  costs  are  much  less.  Major  concerns  are  likely  to  be 
logistics  and  safety  problems  arising  from  the  nitric  acid  component. 

The  specific  blend  studied  In  this  work  was  rather  easy  to  handle 
because  of  the  Intermediate  concentration  of  the  Baume  42  nitric  acid. 
Because  this  blend  seems  to  offer  a  superior  combination  of  cost,  per¬ 
formance,  sensitivity,  and  handling  characteristics,  we  endorse  It  for 
consideration  by  DNA  as  the  baseline  material  in  the  NA/NP  family. 

Recommendation  2:  We  recommend  that  NPN  be  evaluated  for  use  In 
applications  calling  for  an  HE  similar  to  ANFO  but  with  the  ability  to 
do  somewhat  more  work  at  higher  pressures  and  with  more  rapid  kinetics. 

NPN  is  a  good  blasting  agent  with  a  combination  of  properties  that 
may  offer  net  advantages  over  ANFO  in  some  applications.  NPN  Is  similar 
In  many  respects  to  ANFO,  but  has  a  higher  detonation  pressure  and 
velocity,  no  doubt  because  of  Its  greater  Initial  density.  It  also 
probably  has  a  shorter  steady-state  reaction  zone  and  a  more  rapid 
shock-to-detonatlon  (SDT)  transition  than  ANFO.  However,  NPN  costs  more 
than  ANFO,  contains  a  fairly  volatile  solvent  (NP),  and  Is  difficult  to 
keep  In  a  uniform  density  state  because  It  compacts  Irreversibly  In 
response  to  local  pressures.  The  Irreversible  compaction  can  be  used  to 
advantage  by  compressing  the  NPN  to  higher  (presumably  uniform)  densi¬ 
ties  than  studied  here,  where  Its  differences  from,  and  potential  advan¬ 
tages  over,  ANFO  are  greatest. 


NPN  may  be  useful  to  DNA  not  only  In  NWE  simulators,  but  also  as  an 
ANFO  model  In  studies  to  define  the  nonideal  AN  reaction  phenomeno¬ 
logy.  Because  NPN  displays  many  of  the  qualitative  detonation  behavior 
properties  of  ANFO  but  on  a  shorter  time  scale,  It  permits  smaller,  more 
cost  effective  experiments. 

Recommendation  3:  We  strongly  recommend  Lagrange  gage  studies  to 
characterize  SDT  In  NA/NP. 

Rather  slow  nonreproduclble  shock  Initiation  of  detonation  was 
observed  In  one  Lagrange  SDT  experiment  In  NA/NP.  In  addition,  the 
detonation  experiments  Indicated  a  surprisingly  long  steady-state  reac¬ 
tion  zone  In  NA/NP.  Both  phenomena  can  affect  the  design  and  perfor¬ 
mance  of  some  NWE  simulators  and  can  be  well  defined  through  a  series  of 
Lagrange  gage  SDT  experiments.  Such  experiments  are  the  most  direct 
means  we  know  to  characterize  the  reaction  zone  as  It  builds  to  the 
steady-state  and  to  obtain  engineering  curves  of  detonation  buildup 
useful  In  simulator  design  and  modeling. 

Recommendation  4:  We  recommend  standard  rate  stick  studies  In  NPN 
to  determine  the  threshold  diameter  for  size-independent  (Infinite 
diameter)  detonation. 

The  present  NPN  results  are  believed  to  be  a  good  representation  of 
Infinite  diameter  detonation  properties,  but  direct  confirmation  Is 
desirable.  The  most  cost-effective  way  to  obtain  this  Information  Is 
through  rate  stick  experiments  varying  charge  diameter.  The  detonation 
buildup  Information  needed  to  evaluate  or  use  NPN  for  specific  simulator 
applications  can  also  be  generated  In  these  experiments.  Later,  for 
specific  applications,  buildup  distances  can  be  determined  In  dedicated 
experiments  simulating  the  actual  Initiation  conditions,  and/or  Lagrange 
gage  SDT  experiments  can  be  performed  If  more  extensive  understanding  of 
the  buildup  process  Is  required. 

Recommendation  5:  We  strongly  recommend  enlarged  Lagrange  gage 
steady  detonation  experiments,  and  computational  development  as  needed, 
on  NPN  (or  ANFO)  to  characterize  AN  reaction  phenomenology  In  these  HEs. 

One  of  the  most  Important  findings  of  this  study  Is  that  AN-rlch 
HEs  such  as  NPN  and  ANFO,  even  In  large  charges,  have  a  complex  unknown 
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AN  reaction  phenomenology  that  differs  markedly  from  that  of  conven¬ 
tional  HEs  and  current  detonation  models*  Since  ANFO  will  probably  be 
used  In  future  large  NWE  simulators,  our  current  Inability  to  closely 
model  Its  behavior  may  be  of  major  significance  to  DNA.  Therefore,  we 
recommend  that  DNA  initiate  experiments  to  unambiguously  determine  the 
AN  reaction  process  In  this  very  useful  class  of  HE. 

Three  types  of  experiments  available  for  this  purpose  are  cylinder 
tests,  aquarium  tests,  and  Lagrange  gage  experiments*  Although  cylinder 
or  aquarium  tests  provide  expansion  data  to  large  relative  volumes  and 
are  therefore  useful,  a  greatly  enlarged  Lagrange  gage  experiment  offers 
two  major  advantages  in  characterizing  the  detonation  processes  In  these 
HEs*  First,  the  Lagrange  gage  experiment  does  not  Involve  Interaction 
of  the  HE  with  a  container,  so  the  unperturbed  peak  and  final  reaction 
zone  states  characteristic  of  the  explosive  Itself,  rather  than  of  Its 
response  to  the  containment  system,  can  be  determined.  Second,  the 
embedded  instrumentation  resolves  rather  than  Integrates  flow  features 
and  provides  the  Information  necessary  to  Identify  the  termination  of 
the  reaction  and  quantify  Its  rate*  Therefore,  such  experiments,  pos¬ 
sibly  with  NPN  targets  as  suggested  In  Recommendation  2,  are  recommended 
to  generate  a  data  base  characterizing  the  AN  reaction  phenomenology  in 
these  HEs* 

After  the  detonation  process  In  these  AN-rlch  HEs  Is  understood,  it 
may  be  necessary  to  develop  appropriate  computational  techniques  to 
model  It*  Current  DNA  hydrocode  calculations  often  use  the  JWL  equa- 
tlon-of-state  to  represent  the  behavior  of  HE  detonation  products  in 
simulator  applications*  This  EOS  was  originally  developed  to  model  the 
behavior  of  military  explosives,  with  negligible  reaction  zones,  in 
metal  pushing  applications  where  the  roost  significant  energy  exchanges 
occur  near  the  CJ  Isentrope  at  relative  volumes  less  than  about  10*  If 
the  JWL  EOS  Is  used  to  calculate  the  performance  of  ANFO  in  NWE  simula¬ 
tors,  It  may  be  applied  to  much  different  conditions:  long  reaction 
zones,  two-stage  reactions,  low  pressures,  and  states  far  from  the  CJ 
Isentrope. 
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The  effects  of  these  extensions  of  present  computational  techniques 
on  the  accuracy  of  the  calculations  are  unknown.  Therefore,  It  Is 
Important  to  determine  the  sensitivity  of  specific  simulator  design  and 
performance  calculations  to  these  effects  by  parmeter  variation  calcula¬ 
tions  and,  If  necessary,  to  Incorporate  more  detailed  reactive  hydro¬ 
dynamics  or  a  modified  EOS  In  the  hydrocode  calculations. 
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